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INTRODUCTION 


The  Fifteenth  DoD  Tri-Sendee  Review  Conference  on  Atmospheric  Transmission  Models 
was  held  at  the  Geophysics  Directorate,  Hanscom  AFB,  Massachusetts  on  2-3  June  1992.  The 
purpose  of  the  meeting  was  to  review  progress  in  the  modeling  of  radiation  propagating  through 
the  earth’s  atmosphere,  identify  deficiencies  in  these  models,  and  make  recommendations  for 
improvements. 

Approximately  120  scientists  and  engineers,  representing  DoD,  other  government  agencies, 
industry,  and  the  academic  community  were  in  attendance.  The  agenda  consisted  of  twenty  six 
papers  distributed  into  sessions  on:  Code  Enhancements,  Atmospheric  Propagation  Models, 
Measurements  and  Models,  and  Aerosols  and  Clouds, 

This  proceedings  volume  summarizes  the  technical  presentations  at  the  conference.  The  main 
part  of  die  report  consists  of  abstracts  and  copies  of  the  viewgraphs  or  slides  and  other  material 
as  provided  by  the  authors  of  the  presentations.  The  Appendix  includes  the  original  call  for 
papers,  the  invitation,  and  a  copy  of  the  final  agenda.  An  author  index  is  at  the  aid  of  the 
proceedings. 

On  a  personal  note,  this  13th  Annual  Conference  was  the  first  not  organized,  run,  and 
attended  by  Francis  X.  Kneizys.  Hie  "default"  organizers  wish  to  acknowledge  the  impact  of 
his  absence.  While  Frank  took  great  care  to  train  his  apprentices,  his  skills  in  organization, 
maintenance  of  rich  historic  friendships  (including  instant  name-recall),  and  constant  flexibility, 
edged  with  an  intellectually  exciting  (though  subtle)  sense  of  humor,  cannot  be  transmitted 
through  any  traditional  approaches.  The  current  organizers,  along  with  past  and  present 
attendees,  hereby  extend  to  Frank  and  his  family  a  heartfelt  THANK  YOU  and  BEST  WISHES 
for  continued  happiness  and  good  health  in  their  retirement. 
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IMPROVEMENT  OF  MODTRAN  AND  FASCOD  TO  FULL  SOLAR  CAPABILITY 
INCLUDING  MULTIPLE  SCATTERING  AND  SPHERICAL  GEOMETRY 

M.  Ych 

Caelum  Research  Corporation,  1 1229  Lockwood  Drive,  Silver  Spring,  MD  20901 

K.  Stamnes,  S.  Tsay 

Geophysical  Institute,  University  of  Alaska,  Fairbanks,  AK  99775-0800 

We  have  started  a  significant  upgrade  of  MODTRAN  and  FASCOD  to  encompass  full  solar 
capability,  including  rigorous  multiple  scattering  and  spherical  geometry.  As  part  of  this  project 
we  are  attempting  to  modularize  and  restructure  die  existing  codes  to  enhance  the  flexibility  arid 
readability  and  thus  facilitate  their  use  as  well  as  die  inclusion  of  future  improvements.  We  plan  to 
unify  die  structure  of  these  codes  to  provide  a  framework  within  which  die  atmospheric  conditions 
and  optical  properties  other  than  gaseous  absorption  are  a  common  feature.  Thus,  the  desired 
spectral  resolution,  which  dictates  the  treatment  of  gaseous  absorption  by  atmospheric  constituents, 
is  die  key  feature  distinguishing  the  codes.  To  solve  the  ensuing  radiative  transfer  problem  we  will 
utilize  existing  state-of-the-art  multiple  scattering  algorithms.  The  end  product  of  this  effort  will  be 
modem,  well-documented  and  tested  software  tools  that  can  easily  and  safely  be  used  by  others  in 
a  variety  of  applications. 
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UPGRADED  LINE-OF-SIGHT  GEOMETRY  PACKAGE  FOR  MODTRAN 

P.K.  Acharya,  D.C.  Robertson,  A.  Berk 

Spectral  Sciences,  Inc.,  99  South  Bedford  Street,  #7,  Burlington,  MA  01803-5169 

The  geometry  package  in  MODTRAN  was  modified  so  that  die  actual  line-of-sight  (LOS) 
parameters  used  for  the  transmittance/radiance  calculations  always  closely  match  the  user's 
requested  input  values.  For  example,  the  case  specified  by  (HI,  ZENITH  ANGLE,  RANGE) 
would  sometimes  result  in  an  output  range  that  differed  from  the  input  by  20%  or  more.  The 
geometry  routines  incurred  roundoff  problems  on  32-bit  machines  due  to  subtraction  of  large 
numbers,  each  of  which  includes  the  square  of  die  earth-radius.  The  new  routines  replace 
numerically  unstable  algebraic  expressions  by  more  stable  identities  and  occasionally  resort  to 
double  precision  arithmetic.  Furthermore,  die  geometry  routines  can  now  deal  with  very  short 
paths  down  to  0.001  km.  This  new  geometry  package  will  be  part  of  the  next  release  of 
MODTRAN. 
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IR  VISIBILITY 
C.N.  Touart 

Hughes  SIX  Corporation,  109  Massachusetts  Avenue,  Lexington,  MA  02173 

In  response  to  denrand  from  field  userc,dre  current  release  of  the  Phillips  Lab  EOTI>A(EO 
Tactical  Decision  Aid)  outputs  values  of  an  8-12  ^tm  "visibility"  that  is  denoted  IRVIS.  It  is 
defined  as  die  MRT  detection  range  for  an  "avenge"  FLIR  viewing  a  target  of  "standard"  size  and 
inherent  contrast  The  rationale  for  this  definition  and  how  IRVIS  behaves  in  various  weather 
conditions  will  be  presented. 
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A  MODTRAN  ALGORITHM  FOR  PREDICTING  DIRECTIONAL 
SKYSHINE  IRRADIANCE 

D.C  Robertson 

Spectral  Sciences,  Inc.,  1 1 1  S.  Bedford  Street,  Burlington,  MA  01803 

J.A.  Krisd 

Stewart  Radiance  Laboratory,  139  Great  Road,  Bedford,  MA  01730 

B.P.  Sandfond 

Phillips  Labcratory/GPOA,  Hanscom  AFB,  MA  01731-5000 

A  set  of  empirical  algorithms  has  been  developed  to  predict  the  directional  dependence  of 
skyshine  inadianoc  calculated  with  MODTRAN.  They  provide  a  spectral/directional  skyshine 
function  that  allows  skyshine  to  be  predicted  for  arbitrary  directions  from  just  a  few  (four) 
MODTRAN  calculations.  The  algorithms  include  the  thermally  emitted  and  multiply-scattered 
skyshine  radiance  components,  which  are  not  uniform  in  azimuth  and  elevation.  This  approach  is 
applicable  to  many  background  and  target  signature  models  which  often  approximate  the  incident 
skyshine  by  a  contrast  Illustrative  calculations  are  presented  to  show  its  impact  on  calculations  of 
die  reflected  skyshine  signature  components. 


r  i 


L  A 


41 


5:  o 


•-H  C_) 
QC  Uj 
O  QC 
CD  *~-s 


q:  o  ^ 

I- —  * — I  1 _ . 

Q  Q 


Q 
U.  g 

§i* 

TO 

Co  ^  p 

_j  u. 
r>  to 

^S;S 

£;jy 

IIhk 
OQ?  3 


31s 

5"  . 

a  pc 
Uj2g 

OfcS 

JWffi 


C" 

Dr  1° 

J-u  iu 

UJ^ 

OQes 

O 

QC  8  a? 
Qsg 

tr»SS 

^gg 

Q  Sj  m 


CO  i 

hs 

uj  £ 
o  oc 


& 

o  CD 
iu  CO 

5  ^ 

6  Uj 

g:  k, 


oS, . 

S£g 


£g£ 

s  ■* 

£S| 

P5 

51 

^C  iu 


42 


2  JUNE  1992 


* 


o 

1 — i 

I — ■ 
CO 

Q 

Q 

QC 

h- 

1 — i 


CD 

cc 

Lu 


a: 

cn 

cn 

k-  m 

y  § 

O  £ 

£  3 

Q_  3 
O 
O  — J 

2 

5  * 
£  5 

*-i  »— 

cc  Q 
O  § 
CD  * 
— I  3g 

^  s 

2C  Q 

'*  ft 

5 


■ — J  Lu 

o  to 
e>  Q 


ft 

i 

i 


i 


a 

5 

Q 


ID 

2 


i 

CD 

ft 

cn 


ca  in 

i 

in 

£ 


S 


.  co 
to  »-« 
a 

*~i  O 
QC  cc 

t— i  “ 

a  o 
cn  cn 


cn  cn 
cc  *~h 

*— i 

l —  Lu 

cn  2 

w5 

r£ 
cn  ^ 

•  cn 
a j 


cn  Uj 

-J  Q 

LU 

O  CO 
Q  3T 


Uj  J — 

52 

K-  K— 

§  Uj 
CD  ^ 
*-n  ID 

cn  cn 
cn 
cc  ^ 


lT1  ^  7 
1 —  2:  o 

O  s  ^ 

2E  CL 


-C 
lu  £ 

^  s 

cn  o 
_  o 
^  in 

i 

>-  — 
CQ  , 

Q  J 

t  ’ 

H  V- 

><  cc 

Qc  oc 
9- 

^  § 

>-5 


2 

cn  *-~ 

ID 

cn  2 

*-i  CD 

*~s 

Lu  »— 

c  ^ 

H  ~J 

31  p 
cn  co 
>-  -o 

cn  cd 


% 

cc 

*— 

o 

cn 


CD 

5 

CO 

5 


cn 

CD 


cn 


Ui 


§ 
0 D 


45 


ZENITH  ANGLE  (dec) 


K- 

_ x 

lu 

CD 

c o 

co 

Lu 

D* 

CO 

cr  ^ 

CD 

CD 

lu 

3d  Lu 

cr 

CD 

lu 

a 

K- 

CD 

CO 

<c 

CD 

Uj 

CD  Lu 

rU 

lu 

-j 

3D 

Lu 

3; 

CD 

CJ 

QC 

-J  *- 

co  2: 

-J 

3D 

*q:  *-h 

1— — 

»~h  ^  CD 

3: 

Cj 

i —  2 

v-  cr 

^  3:  cd 

*-h  Lu 

CD 

Uj  K*  H 

cr  3: 

> — < 

CO  CD  *■— 

o  K- 

cr  co 

CO 

CO  ^  <C 

CD  CD 

1—  lu 

Lu  U.J  -J 

-j 

CD  — 1 

CD 

_j  3D 

cr 

CD  CD 

h— 

lu  lu  CJ 

CD 

5:  2L 

ISj  — J 

CD  Lu 

•—1 

►i  ^ 

21 

CD 

QC  ^  Cj 

►-H  Co 

h—  1— 

lu 

lu 

~J  lu 

1 — . 

i—  3:  'z. 

^  • — 

lu  Lu 

o  u  *«c 

CJ  5 

cr  cr 

ac 

^  ^  cr 

co  5: 

Uj 

CD 

cr  Lu  •— 

•— 1 

cl  Ll 

LU 

^  CD 

LU  • — • 

5:  iu 

1 — 

X  >—  O 

co  co 

CD  ^ 

2! 

CJ  ^  2£ 

3D  lu 

Cj  CD 

1— 1 

• 

• 

• 

• 

47 

ANGLE  (deg) 


ZENITH  ANGLE  (deg) 


WAVENUMBER  (cm-1) 


O' I  S  O  0  0  S  O-  O  T- 

(0 ‘ I  -  3ie3/dJS^ui)  OllVti 


52 


A 


SKYSHINE  -  NEAR  HORIZON 


o 

o 

o 

cn 


o 

o 

o 

r-» 


o 

o 

o 

in 


o 

o 

o 

cn 


o 
o 
o 

;0t  0OI  ;.0T  g.OT  e.0T  „_0I  g.Ol 


(T_U13/US/2U13/S^^GMTt)  33NVI0VH 


53 


WAVENUMBER  (cm-1) 


r  ^ 


CO 

h- 

Lu 

co 

co 

QC 

h- 

■«* 

Ul 

CC 

K- 

co 

CO 

O 

CC 

^  »— 

U 

P 

■*c 

3:  »-h 

Uj 

k-  cn 

co 

• — «  ►— < 

UJ 

-u 

p 

Ll 

oc  a 

QC 

co 

Lu 

lu 

o 

CO  CO 

30 

co 

CO 

co 

CO 

h— 

Ul 

CO 

CD 

CO 

QC 

CO 

5: 

111 

CD 

CO 

CO 

30 

Q; 

UJ  qc 

K-l 

3 

CO 

QC 

co 

H 

V— I  J — 

CO 

5 

Lu 

a:  co 

QC 

U. 

Cj 

co  co 

»-« 

CO 

hi 

CO 

H 

rr 

Co 

■  "i 

C A 

£ 

^  a: 

K- 

w  # 

CO 

Lu 

co  »- 

U 

CO 

cc 

X 

CO 

^C  Co 

►» » 

H 

CO 

co  co 

QC  CO 

K» 

L. 

H 

r«k 

*— 1 

CO  QC 

Li  i 

=u 

U 

3:  co 

QC  »"-i 

-J 

CD 

30 

s 

1— -  K~ 

30 

QC 

Co 

CO 

-o 

Q 

u.  *-i 

Uj 

-J 

K* 

o 

Uj  Co 

co 

co 

1 — ■  co 

2 

U0 

Li  ■ 

UJ 

1 

^  Ul 

■**  , 

Cj 

co 

1 

2 

CO  K- 

-J  CO 

CO 

l-H 

H 

K-,  ^ 

30  ^ 

►-H 

g 

QC 

CO  cc 

CO  *-i 

QC 

ac  co 

— I  3C 

& 

,2 

2: 

Uj  Lu 

<c  Uj 

fS. 

UJ 

LU 

Z>  h- 

co  ^ 

Lu 

_ 1 

CO 

2! 

CO 

^  J 

co 

co 

CO 

ct 

cc 

Uj  * — 

1  ^ 

Q  > 

lt! 

•  • 

lt! 

1—  Uj 

CO  Uj 

•  • 

£ 

co 

£ 

»— i  Uj 

CO  QC 

Co 

(/) 

K* 

CD 

2  co 

LU 

QC 

►-H 

3C  U. 

•-H 

QC 

CO 

Ll. 

J — 

22  ^ 

CO  *-1 

CO 

CL 

a: 

CD  CO 

CO 

co 

• 

• 

• 

• 

54 


L.  A 


TARGETS  AND  SCENES  MODEL 
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OBSERVER  ASPECT  ANGLE  (deg) 
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STATUS  OF  PL/GP  HIGH  RESOLUTION 
RADIANCE-TRANSMITTANCE  MODEL:  FASCODE 

EX.  Kneizys,  G.P.  Anderson,  J.H.  Chetwynd,  L.W.  Abreu,  M.L.  Hoke 
Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  AFB,  MA  01731 

S.A.  Clough,  RJD.  Worsham 

Atmospheric  and  Environmental  Research,  Inc.,  840  Memorial  Drive,  Cambridge,  MA  02139 

EJ>.  Shettle 

Naval  Research  Laboratory,  Code  6322,  Washington,  DC  20375 

FASCODE,  a  line-by-line  atmospheric  radiance-transmittance  code,  has  been  released  in  a 
preliminary  version,  FASOOD3P  (Much  1992).  As  with  previous  versions,  FASCOD3P  is 
applicable  to  spectral  regions  from  the  microwave  to  the  middle  ultraviolet,  employing  standard 
spectroscopic  parameter  supplied  from  external  line  atlases.  FASCOD3P  is  the  only  DoD  code 
currently  implementing  all  HTTRAN92  data,  including  the  temperature-dependent  CFC  absorption 
cross  sections.  Other  important  upgrades  include:  weighted  radiance  algorithms,  multiple  line 
coupling  options,  efficient  weighting  functions,  LTE/NLTE  line  shape  compatibility,  and  enhanced 
output  options. 
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HITRAN  1992:  M3  (MULTI-MEDIA  MILIEU) 

L.S.  Rothman 

Geophysics  Directorate,  Phillips  Laboratory/Optical  Environment  Division  Simulation  Branch 

(OPS),  Hanscom  AFB,  MA  01731-5000 

A  new  edition  of  die  spectroscopic  molecular  database,  HITRAN,  became  available  in 
March  1992.  This  current  edition  contains  over  70  Megabytes  of  high  resolution  data  of  transitions 
for  31  species  and  their  atmospherically- significant  isotopic  variants.  In  addition,  there  is  a  file  of 
cross-sections  for  heavy  molecular  species,  with  bands  at  several  representative  temperatures  that 
can  be  used  for  atmospheric  retrievals  with  some  codes  such  as  FASCODE3P.  HITRAN  is  now 
available  on  alternative  media:  floppy  diskettes  and  optical  CD-ROM.  Unlike  the  mainframe- 
standby  nine-track  magnetic  tape,  these  media  permit  rapid  access  and  new  user  features.  This  talk 
will  summarize  some  of  the  major  updates,  improvements,  and  modifications  to  HITRAN  as  well 
as  the  new  PC-oriented  media  options. 
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HITRAN  1992: 

M3  (Multi-Media  Milieu) 

Laurence  S.  Rothman 
Geophysics  Directorate 
Optical  Environment  Division 
Simulation  Branch  (OPS) 
Hanscom  AFB,  MA  01731-5000  USA 


Abstract 

A  new  edition  of  the  spectroscopic  molecular  database,  HITRAN, 
became  available  in  March  1992.  This  current  edition  contains  over  70 
Megabytes  of  high  resolution  data  of  transitions  for  31  species  and  their 
atmospherically-significant  isotopic  variants.  In  addition,  there  is  the  file 
of  cross-sections  for  heavy  molecular  species,  with  bands  at  several 
representative  temperatures  that  can  be  used  for  atmospheric  retrievals 
with  some  codes  such  as  FASCODE3P.  HITRAN  is  now  available  on 
alternative  media:  floppy  diskettes  and  optical  CD-ROM.  Unlike  the 
mainframe-standby  nine-track  magnetic  tape,  these  media  permit  rapid 
access  and  new  user  features.  This  talk  will  summarize  some  of  the 
major  updates,  improvements,  and  modifications  to  HITRAN  as  well  as 
the  new  PC-oriented  media  options. 


Details  (of  both  1991  and  1992  issues)  to  appear  in 
Journal  of  Quantitative  Spectroscopy  and  Radiative 
Transfer,  Volume  48,  1992 

Media  options: 

Magnetic  Tape  (9-track,  6250  PI) 

CD-ROM 

Floppy  diskettes 

HAWKS  - 

Faster  publication 

Supplemental  files  (CFC's;  High-vib  bands) 
Next  HITRAN  Conference  -*  June  1 993 
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THE  HITRAN  MOLECULAR  DATABASE: 
EDITIONS  OF  1991  AND  1992 


L.S.  ROTHMAN1 

Phillips  Laboratory,  Geophysics  Directorate 
Hanscom  AFB,  MA  01731,  U.S.A. 

R.R.  GAMACHE 

University  of  Massachusetts  Lowell,  Center  for  Atmospheric  Research 
450  Aiken  Street,  Lowell,  MA  01854,  U.S.A. 

R.  H.  TIPPING 

Department  of  Physics  &  Astronomy,  University  of  Alabama 
Tuscaloosa,  AL  35487,  U.S.A. 

C.P.  RINSLAND  and  M.A.H.  SMITH 
Chemistry  &  Dynamics  Branch,  NASA  Langley  Research  Center 
Hampton,  VA  23665,  U.S.A. 

D. CHRIS  BENNER  and  V.MALATHY  DEVI 
Department  of  Physics 
College  of  William  and  Mary 
Williamsburg,  VA  23187,  U.S.A. 

J.-M.  FLAUD.  C.  CAMY-PEYRET,  and  A.  PERRIN 
LPMA,  CNRS,  University  P.  et  M.  Curie 
4  Place  Jussieu,  75252  Paris,  France 

A.  GOLDMAN 

Department  of  Physics,  University  of  Denver 
Denver,  CO  80208,  U.S.A. 

S.T.  MASSIE 

National  Center  for  Atmospheric  Research 
P.O.  Box  3000,  Boulder,  CO  80307,  U.S.A. 

L.R.  BROWN  and  R.A.  TOTH 
Jet  Propulsion  Laboratory 

4800  Oak  Grove  Drive,  Pasadena,  CA  91109,  U.S.A. 


*To  whom  all  correspondence  should  be  addressed. 

(to  appear  in  JQSRT,  vol.  48,  1992) 
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New  or  Modified  Parameters  by  Edition 


1991 

1992 

h2o 

★ 

CN 

o 

u 

★ 

★ 

03 

★ 

CO 

★ 

ch4 

★ 

★ 

NO 

★ 

S02 

★ 

N02 

★ 

HN03 

★ 

HF 

★ 

★ 

HCI 

★ 

HBr 

★ 

HI 

★ 

OCS 

★ 

n2 

★ 

h2o2 

★ 

★ 

c2h2 

★ 

c2h6 

★ 

cof2 

★ 

★ 

sf6 

★ 

h2s 

★ 
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Molecule 

Spectral  Range 

(cm'* 1) 

f  Number  of 
s  points  per  T 


CCI3F  (CFC-11) 


830-860 


1060-1107 


CCIjFj  (CFC-12) 

867-937 

1080-1177 

CCIFj  (CFC-13) 

765-805 

1065-1140 

1170-1235 

C2a3F3  (CFC-11 3) 


C2CI2F4(CFC-114) 


C2CIFS  (CFC-11 5) 


780.5-995 


1005.5-1232 


815-860 


870-960 


1030-1067 


1095-1285 


955-1015 


1110-1145 


1167-1259 


n2o5 

555.4-599.8 

720.3-764.7 

1210.1-1274.8 

1680.2-1764.6 

„  430 

454 


3  034 


12  808 


CI0N02 

740-840 

1240-1340 

1680-1790 

770-830 


786-806 


1255-1290 


10  371 


5  476 


CHCI2  (CFC-21) 

785-840 

1 

5  020 

1  826 


2  359 


11  798 


7  301 


CF4  (CFC-14)  1 


CHCIF2  (CFC-22)  1  780-1335 


1270-1350 


1  Omitted  from  HITRAN'9l,  added  to  HtTRAN'92. 

1  Removed  from  HITRAN.^tt  (already  on  high  resolution  portion). 
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1992  CD-ROM 


HITRAN  92  CD-ROM 
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Future  Version  1993 
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RAN-PC 
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Smaller  Random 
Access  Data  File 
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per  day,  7  days  a  weeK  beginning  i  June 

Problem  reports  should  be  directed  to  the  System  Operators 
Dale  Sinclair,  377-3694  Bob  Hawkins,  377-3614 
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STATUS  OF  PL/GP  MODERATE  RESOLUTION  RADIANCE/ 
TRANSMITTANCE  MODELS:  MODTRAN  AND  DERIVATIVES 

FX  Kneizys,  L.W.  Abrcu,  G.P.  Anderson,  J.H.  Chetwynd 
Geophysic  Directorate,  Phillips  Laboratory/GPOS,  Han  scorn  AFB,  MA  01731 

MODTRAN  (public  release:  November  1990)  is  a  2  parameter  (pressure  and  temperature) 
band  model  code  with  moderate  spectral  resolution  (2cm  **  full  width-half  maximum).  The 
MODTRAN  band  model  parameters  ere  calculated  by  directly  utilizing  the  HTTRAN-86  data 
base.  The  resultant  code  is  fully  compatible  with  LOWTRAN  7  capabilities  and  has  been 
successfully  validated  against  FAS  CODE. 

Because  of  its  relative  accuracy  and  ease  of  use,  MODTRAN  is  undergoing  an  applications 
revolution.  For  instance,  it  serves  as  a  primary  foundation  for 

AURIC-E.  die  Atmospheric  Ultraviolet  Radiance  Integrated  Code  (Extension),  describing  solar 
scattering  and  transmittance  to  0.15  microns. 

SAM.  SHARC  and  MODTRAN,  a  newly  envisioned  rapid  combined  algorithm  far  surface  to 
space  NLTE  calculations. 

MQSART.  Moderate  Spectral  Atmospheric  Radiance  and  Transmittance  Code,  a  combination  of 
the  MODTRAN  band  model  and  the  APART  scene  and  cloud  formulations. 

FAUST.  Fast  Atmospheric  Unified  Synthesis  Technique,  a  pragmatic  limb  viewing  algorithm, 
combining  MODTRAN  and  SHARC  path  segments. 
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1992  Annual  Review  Conference  on 
Atmospheric  Transmission  Models 


STATUS  OF  PL/GP  MODERATE  RESOLUTION 
RADIANCE/TRANSMITTANCE  MODELS: 

MODTRAN  AND  DERIVATIVES 


F.X.  Kneizys,  L.W.  Abreu,  G.P.  Anderson,  and  J.H.  Chetwynd 
Geophysics  Directorate,  PtyGPOS,  Hanscom  AFB 


MODTRAN  (public  release:  November  1990)  is  a  2  parameter 
(pressure  and  temperature)  band  model  code  with  moderate 
spectral  resolution  (2ca-l  full  width-half  maximum) .  The 
MODTRAN  band  model  parameters  were  calculated  by  directly 
utilizing  the  HITRAN-86  data  base.  The  resultant  code  is  fully 
compatible  with  LOWTRAN  7  capabilities  and  has  been  successfully 
validated  against  FASCODE. 


Because  of  its  relative  accuracy  and  ease  of  use,  MODTRAN  is 
undergoing  an  applications  revolution.  For  instance,  it  serves 
as  a  primary  foundation  for: 

AURIC-E .  the  Atmospheric  Ultraviolet  Radiance  Integrated 
Code  (Extension) ,  describing  solar  scattering  and 
transmittance  to  0.15  microns. 

SAM.  SHARC  and  MODTRAN,  a  newly  envisioned  rapid  combined 
algorithm  for  surface  to  space  1LTE  calculations. 

MOSART.  Moderate  Spectral  Atmospheric  Radiance  and 
Transmittance  Code,  a  combination  of  the  MODTRAN  band 
model  and  the  PRA/ APART  scene  and  cloud  formulations. 
FAUST .  Fast  Atmospheric  Unified  Synthesis  Technique,  a 
pragmatic  limb  viewing  algorithm,  combining  MODTRAN  and 
SHARC  path  segments. 
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GPOS  RADIATIVE  TRANSFER  CODES: 


LOWTRAN 
MODTRAN 
FAS CODE 


Low  Spectral  Resolution  Transmission  and  Radiance  (LOWTRAN) 

Code:  , 

LOWTRAN  (L0WTRAN7,  1989)  Is  a  low  spectral  resolution  (20cm”1) 
code  for  predicting  the  atmospheric  transmittance  and  background 
radiance  for  line-of-sight  paths  at  altitudes  below  50  km  over 
the  spectral  range  from  0  to  50,000  cm”1.  The  radiance 
includes,  in  addition  to  emission  from  13  atmospheric  gases, 
contributions  from  single  scattered  solar  and  lunar  radiance  and 
multiple  scattered  solar  and  thermal  radiation.  Direct  solar 
irradiance  is  also  an  option.  A  single  parameter  band  model  is 
used  for  molecular  line  absorption  (based  on  spectrally  degraded 
"exact"  line-by-line  calculations  using  FASCOD2  (see  below) ,  and 
the  HITRAN86  data  base) .  Molecular  continuum-type  absorption 
spectra  (for  H20,  C02,  02,  N2,  and  03)  have  been  developed 
separately  and  subsequently  incorportated  into  LOWTRAN. 

Molecular  scattering,  aerosol  and  hydrometeor  absorption  and 
scattering,  and  surface  emittance/reflectance  properties  are 
also  included,  as  well  as  a  self-contained  spherical  viewing 
geometry  algorithm  that  includes  the  effects  of  refraction  and 
earth  curvature.  Representative  atmospheric  gas,  aerosol,  cloud 
and  rain  models  are  provided  with  options  to  replace  them  with 
user-defined  values.  This  fast-running  code  has  long  been  the 
standard  for  broad-band  electro/optical  systems  design  in  the 
ultraviolet  through  infrared  frequencies.  With  the  extension  to 
the  microwave,  the  code  can  also  be  employed  as  a  primitive 
energy  balance  (climate)  code.  Future  options  will  center  on 
maintaining  parallelism  with  MODTRAN;  see  below. 


Moderate  Resolution  Uanalaalan  and  Radiance  fiate:  (modtran) 

MODTRAN  (MODTRAN,  1990)  is  a  moderate  resolution  (2  cm”1)  code 
for  predicting  the  atmospheric  transmittance  and  background 
radiance,  appropriate  for  altitude  regions  not  effected  by  non¬ 
local  thermodynamic  equilibrium  (non-LTE)  signatures,  generally 
below  60km.  A  two-parameter  (temperature  and  pressure)  band 
model  is  used  for  molecular  line  absorption,  based  directly  on 
the  HITRAN86  data  base,  employing  FAS CODE  exact  calculations  as 
validation  only.  MODTRAN  contains  all  the  options  and  spectral 
coverage  (0  *o  50,000  cm”1)  of  LOWTRAN  7,  including  default 
gases,  aerosol,  cloud  and  rain  profiles,  surface  properties, 
continue  contributions,  viewing  geometry,  and  solar/lunar  and/or 
thermal  direct  or  multiply  scattered  radiance.  User-defined 
input  can  replace  default  options.  Because  of  its  improved 
spectral  resolution  and  accuracy,  MODTRAN  is  recommended  as  the 
appropriate  replacement  for  LOWTRAN  7  for  systems  evaluation. 
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including  transmission  of  aircraft  and  pluae  signatures. 

Suitable  scanning  and  filter  options  are  available  to  natch  E/0 
systems  resolution.  Note  that  KODTRAN  explicitly  contains 
LOWTRAN  7  as  an  operational  node  to  facilitate  comparisons. 
Future  options  will  include  addition  of:  enhanced  UV  capability 
(AURIC),  man-made  contaminant  signatures  (CPC's) ,  suitable 
modifications  for  climate  simulations,  and  improved  surface  and 
cloud  properties  (MOSART) . 


Fflgt  Afagcharig  Signature  Code  (fas code) 

FAS CODE,  most  recently  released  as  FASC0D3P,  a  preliminary 
version  of  FASC0D3  (1992) ,  is  a  state-of-the-art  line-by-line 
transmittance  and  radiance  simulation  code.  For  atmospheric 
cases,  it  has  all  the  flexibility  and  default  options  (UV  to  mm- 
wave,  clouds,  rain,  surface  properties,  etc.)  of  LOHTRAN  and 
KODTRAN,  with  the  exception  of  solar/lunar  irradiancea.  The 
altitude  range  is  essentially  unlimited  (0-300km)  but  default 
constituent  profiles  are  currently  specified  only  to  120km. 
Because  it  explicitly  calculates  the  exact  line  shape  for  the 
set  of  species  specified  by  the  HI TRAN  spectroscopic  parameters 
adjusted  for  local  temperature  and  pressure,  FASCODE  can 
accurately  follow  arbitrarily  configured  lines  of  sight, 
including  those  that  cross  from  Doppler  to  Voigt  regimes  and/or 
intersect  the  surface.  The  code  also  has  a  non-LTE  capability, 
incorporating  contributions  from  non-thermal  emission  sources, 
given  appropriate  populations  (as  precalculated  by  SHARC, 
below) .  FASCOD3P  is  the  only  DoD  code  that  has  fully 
implemented  the  fluorocarbon  cross  sections  in  conjunction  with 
local  heating/cooling  rate  estimates,  enhancing  its  role  as 
validation  standard,  climate  code,  and  "forward"  radiance 
algorithm  for  remote  sensing.  Operational  speed  of  FASCODE  can 
be  prohibitive  for  real  time  analyses,  but  the  accuracy  is 
essential  for  evaluating  more  pragmatic  algorithms.  Future 
endeavors  include  additions  of:  full  vectorization,  improved 
line  shape  for  remote  sensing,  solar  capability  with  new 
multiple  scattering  algorithm,  and  optimized  feedback  to 
radiative/photochemical/dynamical  programs. 


Mvdvrate  spectral  Af ggphsrlc  Raflianss  and  Uanaltfcancft*  mosart 

MOSART  is  a  "spatially  extended"  version  of  M0DTRAN,  providing  much 
enhanced  descriptions  of  constituent,  aerosol,  cloud  and  surface 
properties.  The  new  data  bases  and  implementation  algorithms  for 
the  complex  scene  descriptors  will  be  derived  from  the  APART  Code 
(Photon  Research  Associates)  but  the  spectroscopic  foundation  will 
remain  KODTRAN.  All  ongoing  KODTRAN  enhancements  will  be 
maintained.  This  new  code  is  expected  to  readily  into  SSGM. 


r  i 


MODTRAK  Authorship 


F.X.  Kneizyp,  L.W.  Abreu,  E.P.  Shettle1, 
J.H.  chetwynd,  G.P.  Anderson,  w.o.  Gallery2, 
J.E.A.  Selby3,  S.A.  Clough2 
(PL/ Geophysics  Directorate) 

and 

A.  Berk,  L.S.  Bernstein,  D.C.  Robertson 
(Spectral  Sciences,  Inc.) 


POC/LOWTRAN  and  MODTRAN 


Gail  P.  Anderson 
PL/GPOS 

Hanscos  AFB,  MA  01731 


617-377-2335 
617-377-2337 
FAX  617-377-8780 


AVAILABILITY 

NATIONAL  CLIMATIC  DATA  CENTER  (NOAA) 
CLIMATE  SERVICE  SECTION 
FEDERAL  BUILDING 
ASHEVILLE,  NC  28801 

704-259-0272 


or 

DoD  Reguests  to  G.P.  Anderson 
as,  above 
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1992  PLANS 
MODTRAN:  Evolution, 


Documentation,  Mutations 


MODTRAN 

NEAR  TERM  - 


MODTRANn  AURIC  MOSART  SAM  etc. 


Geometry/ Precis ion  *SSI 
03  9.6  micron  band  *SSI 
1992  HITRAN  database  *SSI 


1R  CFC's 


? 


uv: 

02  S/R 
CFC'S 
S02,  N02 

Vis: 

03  Chappuis  ? 

Multiple  Scattering  Caelum 

Very  Low  sun  Solution  AER? 


Aerodyne 

? 

? 


PHOTOCHEMISTRY-related  - 

Heating/Cooling  Rates 
FSC3  Validation 

Photodissociation  Rates  ?  ? 

Climate  Change  Scenarios  ? 


INVERSION  ALGORITHM  DEVELOPMENT  - 


Inversion  Algorithm: 

Feasibility  ? 

FSC3  Validation  ? 


SENTRAN  Extensions  SPARTA 

A-matrix  or  vgt.fun.  ? 
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PATH  DESCRIPTORS 


Cl imatologies : 

Cloud,  Aerosol  ? 

CPC's  NRL 

I sage  Properties:  APART 
Surface 
Clouds 
Backgrounds 


SSI/PRA 

n  j  m 

*  /  * 


OTHER  - 


Compatibility 


SHARC 

SSGM 

SMAART 

FAUST 


Cooperative  Efforts 
NASA 


Model  Comparisons 
Ozone  Depletion 
ITRA/ICRCCM 


SSI 

PRA/? 

MRC/? 

MRC/? 


Harvard 

APL 

DOE/ WHO 


L  J 


81 


r  i 


THE  ATMOSPHERIC  ULTRAVIOLET  RADIANCE  INTEGRATED  CODE 

(AURIC) 

GC.  Betchley,  J.A.  Con  ant 

Aerodyne  Research,  Inc.,  45  Manning  Road,  Billerica,  MA  01821 
GJP.  Anderson,  L.A.  Hall 

Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  AFB,  MA  01731 

This  project  is  a  multi-year  effort  to  develop  new  capabilities  for  simulating  background 
radiance  and  transmission  throughout  die  atmosphere,  as  viewed  by  ultraviolet  sensors.  The  recent 
focus  has  been  to  develop  an  initial  version,  built  upon  the  existing  MODTRAN  model.  This 
version  extends  the  altitude  range  from  120km  up  to  1000km,  and  the  wavelength  range  from 
200nm  down  to  120nm. 

The  routines  in  AURIC  were  generalized  to  handle  more  atmospheric  layers  than  with 
MODTRAN.  The  current  dimensioning  is  for  100  layers,  but  can  be  changed  by  editing  one  file 
and  then  re-compiling  the  software.  The  profile  data  in  MODTRAN  were  then  extended  to 
1000km  altitude,  using  data  from  the  SHARC  code  from  120  km  to  300  km  and  from  the  US 
1976  Standard  Atmosphere  from  300  km  to  1000  km.  The  Schumann-Runge  O2  continuum  and 
band  system  model  was  upgraded  using  die  model  developed  by  K.  Minschwaner  of  Harvard 
University.  Coefficients  for  a  set  of  quadratic  expressions  in  temperature  were  derived  to  generate 
appropriate  cross-sections  at  a  resolution  of  1 .0  cm'1.  Recent  SUSIM  solar  data  were  added  in  the 
120  to  400  nm  wavelength  region. 
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Wavenumber  (cm*1)  in  thousands 
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SYSTEMS 


SAM:  SHARC  AND  MODTRAN  MERGED 


A.  Berk,  D.C.  Robertson,  L.S.  Bernstein 
Spectral  Sciences,  Inc.,  1 1 1  S.  Bedford  Street,  Burlington,  MA  01803 

R.  Sharma 

Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  AFB,  MA  01731 

An  effort  to  merge  SHARC  (50-300km)  and  MODTRAN  (0-1 20km)  into  a  single, 
seamless,  and  fully  correlated  radiance  code  dubbed  SAM  has  begun.  SAM  will  employ  a  single 
curve-of-growth,  include  all  MODTRAN  features,  and  smoothly  transitions  into  NLTE  conditions 
as  they  arise  above  SO  km.  The  MODTRAN  aerosol  models  will  be  incorporated  into  the  NLTE 
radiance  calculation,  and  an  option  to  use  multiple  profiles  far  a  single  calculation  will  be  made 
available.  An  equivalent  width  radiative  transport  algorithm  has  been  developed  which  is 
applicable  throughout  the  0- 300km  altitude  regime.  Preliminary  results  will  be  presented. 
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REQUIRE  AGREEMENT  WITH  SHARC  ABOVE  50  KM 
AND  M0DTRAN/FASC0D  BELOH  50  KM 


NIGHTTIME  RADIANCE 


WAVELENGTH  (jjtm) 
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SCA TTERING 

-  MQDTRAN  MULTIPLE  SCATTERING  BELOH  50  *'• 

-  SINGLE  SCA  TTERING  ABOVE  50  KM 
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EMPLOYS  REFRACTIVE  PATH  GEOMETRY 
ALL  Oh/S  3-D  ATMOSPHERE 

RADIANCE  PREDICTIONS 


RADIATION  TRANSPORT 
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AVOIDS  LINE  STRENGTH  DISTRIBUTION  ASSUMPTION 
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AN  INITIAL  VERSION  OF  SAM  IS  SCHEDULED  FOR  DELIVERY 
TO  PHILLIPS  LABORATORY  IN  EARLY  1993. 


SMART  MODEL  FOR  ATMOSPHERIC  AND  ASTRONOMICAL  RADIANCE  AND 

TRANSMITTANCE  (SMAART) 


S.B.  Downer,  J.P.  Kennealy,  P.C.F.  Ip 
Mission  Research  Corporation,  One  Tara  Blvd,  Suite  302,  Nashua,  NH  03062 

F.O.  Clark 

Geophysics  Directorate,  Phillips  Laboratory,  Han  scorn  AFB,  MA  01731 

Dramatic  changes  are  occurring  in  defense,  requiring  easily  accessible  background  radiance 
descriptors  for  downward  looking,  TMD,  and  GPALS  scenarios.  A  knowledge-based  solution  is 
offered  which  fills  these  needs  by  combining  die  equilibrium  lower  atmosphere  model  and  non¬ 
equilibrium  upper  atmosphere  model  into  a  unified  seamless  single  answer,  SMAART. 

SMAART  includes  an  easily  accessible  guide  for  users  in  both  the  proper  selections  of  the 
codes  and  their  execution.  It  implements  a  single  interface  to  models  with  rapid  response  validated 
data  bases  to  permit  rapid  exploration  of  parameter  space,  based  on  DoD/SDIO  Standard  codes. 
SMAART  will  also  contain  DoD  Celestial  Models,  and  runtime  capabilities  for  MODTRAN, 
SHARC  and  other  PL  Atmospheric  Models.  Accompanying  this  intelligent  front  end  is  an 
extensive  set  of  data  analysis  tools  that  will  allow  users  to  readily  examine  and  compare  results  for 
any  of  die  SMAART  modeling  codes.  Present  and  future  capabilities  are  presented  along  with  a 
demonstration  of  the  first  prototype,  SMAART  Version  0.5. 
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Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  AFB,  MA  01731 
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Graphically  Interactive  Area-oMnterest  Definition 
Provides  Intensity  Information  at  any  User-Specified  Point 
Provides  Integrated  Intensities  Over  User-Defined  Bands 
Interactive  Image  Display 
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CD  ROM  Release 

More  Extensive  Data  Analysis  Tools 


ONCORE:  THE  ONTAR  IMPLEMENTATION  OF  LOWTRAN  7 


J.  Schroeder 

Ontar  Corporation,  129  University  Road,  Brookline,  MA  02146 

Oncore  [tm]  is  a  user-friendly  implementation  of  the  Phillips  Laboratory/  Geophysics 
Directorate's  LOWTRAN  7  atmospheric  transmission  and  radiance  model  The  software  will 
operate  on  IBM  (or  compatible)  personal  computers  in  die  MS-DOS  environment,  or  the 
workstation  UNIX  environment 

The  oncote  package  includes  software  for  selecting  inputs  -  with  on-line  help,  viewing 
results  on  die  computer  screen,  and  producing  hard  copy  outputs.  Additional  features  include 
utilities  for  exporting  compatible  files  for  many  word  processing  and  desk  top  publishing 
packages. 

This  presentation  will  demonstrate  the  oncost  software,  and  discuss  sample  calculations. 


ONCORE  [c]:  An  implementation  of  the  PL/GP's 
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Brookline,  MA  02146-4532 

Tel:  617-739-6607  FAX:  617-277-2374 


PCTRAN  7  [c]  *  Version  2a 
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Setup  Program 


ONTAR  Corporation,  129  University  Road,  Brookline,  MA02 
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MICROWAVE  REFRACTIVE  PATH  CALCULATIONS  FOR  THE  PHILLIPS 
LABORATORY  (GL)  RADTRAN  ATTENUATION/TRANSMITTAN 
CE/BRIGHTNESS  TEMPERATURE  COMPUTER  CODE 

R.G.  Isaacs,  W.O.  Gallery,  R.D.  Worsham 

Atmospheric  and  Environmental  Research,  Inc.,  840  Memorial  Drive,  Cambridge,  MA  02139 

VJ.  Falcone 

Geophysics  Directorate,  Phillips  Laboratory,  Hanscotn  AFB,  MA  01731 

The  RADTRAN  computer  code  was  developed  by  the  Air  Force  Geophysics  Laboratory  (GL)  to 
provide  atmospheric  attenuation  and  brightness  temperature  calculations  fbrtypical  atmospheric 
paths  over  the  frequency  range  from  1  to  300  GHz  (Falcone  et  al.,  1982).  The  atmospheric 
attenuation  submodels  of  the  clear  atmosphere,  fog,  cloud,  and  rain  used  in  RADTRAN  have  been 
thoroughly  documented  (Falcone  et  aL,  1979).  The  enhanced  RADTRAN  is  described  in  Isaacs  et 
al  (1989).  In  addition  to  calculation  of  attenuation,  transmittance,  and  brightness  temperature,  this 
version  treats  surface  emissivity  (Isaacs  et  al.,  1988),  scattering  properties  of  precipitation  (Isaacs 
et  al,  1987),  and  multiple  scattering  (Jin  and  Isaacs,  1987). 

In  this  paper  we  describe  an  additional  enhancement  to  include  the  curved  earth  and 
refraction  in  the  determination  of  path  properties.  The  atmospheric  ray  trace  algorithm  is  based  on 
the  PL  Optical  Physics  Division  FSCASTM  program  (Gallery  et  al.,  1983)  which  calculates  the 
integrated  absort>er  amounts  for  FASCODE.  This  algorithm  has  been  adapted  to  the  microwave 
region  and  is  used  to  calculate  the  path  absorber  amounts  for  the  attenuation  and  transmittance 
modes  of  RADTRAN.  Applications  to  high  incidence  angle  viewing  will  be  discussed. 
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The  RADTRAN  Computer  Code 
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Falcone  et  al,  1982  { Proc .  Soc.  Photo.  Opt  tnstrum.  Eng. 
Isaacs  et  al,  1989  { GL-TR-89-0188 } 
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Spherical  Refractive  Geometry,  same  as  in  LOWTRAN  and 
FASCODE  {Gallery  et  al.,  1983,  AFGL-TR-83-0065} 
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Geometry  of  the  refraction  calculation 
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Refraction  Effects,  Sensor  at  Ground 


Refraction  Effects,  Sensor  at  Ground  (con t) 


Zenith  Angle  (deg)  Zenith  Angle  (deg) 
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ATMOSPHERIC  TRANSMISSION  MODELS  FOR  PASSIVE  MILLIMETER 
WAVE  TARGET  SIGNATURE  SIMULATION 

W.T.  Kreiss 

Horizons  Technology,  Inc.,  9390  Ruffin  Road,  San  Diego,  CA  92123-1826 

R.S.  Dummer 

R.S.  Dummer  and  Associates,  Inc.,  San  Diego,  CA  92131 

Passive  Millimeter  Wave  (MMW)  sensors  have  been  seen  as  an  adjunct  to  VISEO  sensors 
for  multi-spectral  surveillance,  target  detection,  and  tracking  for  several  decades  due  to  their 
superior  obscuram  penetrating  capability.  The  "antenna  proMem"  has  been  the  main  deterrent  in 
their  successful  application.  Current  work  on  sensor  front-end  design,  and  on  improved  signal 
processing  algorithms,  are  keeping  die  promise  of  combined  MMW/VISEO  systems  alive. 
However,  a  number  erf  problems  associated  with  MMW  atmospheric  transmission  and  radiance 
modeling  still  remain  to  be  resolved.  For  this  reason  several  existing  MMW  tnutsmissfoo/Sradiance 
models,  including  LOWTRAN  7,  have  been  examined  to  assess  their  capabilities  for  addressing 
these  problems.  This  presentation  will  describe  targeting  scenarios,  problems  associated  with 
atmospheric  structure  and  composition,  and  the  ability  of  several  MMW  transmissionfradiancc 
models  to  meet  the  requirements  of  aerial  targeting  applications. 
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APPARENT  BRIGHTNESS  TEMPERATURE  OF  SUBTROPICAL  OCEAN  SURFACE 

FROM  10  KM 
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BRIGHTNESS  TEMPERATURE  OVER  SUBTROPICAL  OCEAN  FROM  10  KM 
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Receiver  Internal  Noise  Temperature  (k) 


Critical  Areas  In  Phenomenology 
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Fixed  Target  Detection 
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Moving  Target  Detection 
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35  GHz  Antenna  Diameter  -  (Time-Bandwidth1)  Tradeoffs  For 
Target  Cross  Sections 


Oiameter  ol  Circular  Aperture  (Meters) 
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95  GHz  Antenna  Diameter  -  (Time-Bandwidth)  Tradeoffs  For 
Target  Cross  Sections 
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EFFECT  OF  CLOUDS  ON  MILLIMETER  WAVE  BRIGHTNESS  TEMPERATURES 


I 


MODEL  =  BACKFIRE 
ALT  =  iOKM 
HEADING  =0.0 

YAW  =  PITCH  =  ROLL  =  0.0  DEG 
EARTH  CENTRAL  ANG  =  0.0  DEG 
TGT  TEMP  =  250  K 
SENSOR  ALT  =  0.0  KM 
CENTRAL  OBS  FREQ  =  35  GHZ 
BANDWIDTH  =  1  GHZ 
ATMO  =  US  STD 
SOLAR  IRRAD  -  0.0  W  M**2 


SURFACE  =  LAND 
SFC  TEMP  =  290  K 
SFC  TYPE  =  VEGETATION 
SFC  FRACTIONAL  COVER  =  0.75 
SFC  DIELECTRIC  CONST  =  10 
SFC  BACKSCATT  COEF  =  02 
CLOUD  =  NONE 
CLOUD  BOTTOM  = 

CLOUD  THICK  = 

CLOUD  LWC  = 


OBSERVER  DISTANCE  -  10  KM 
MEAN  ATMO  XMISSION  =  0.9350 
TGT  PROJ  AREA  -  290.9  SQ  M 
BACKGND  TEMP -21.1  K 
MEAN  TGT  TEMP  =  249.7  K 
MEAN  APPARENT  TEMP  =  213-7  K 
MEAN  INHERENT  TEMP  «  228.6  K 
TGT  APPARENT  CONTRAST  =  10.13 
TGT  INHERENT  CONTRAST  =  10.83 
APPARENT  /  INHERENT  RATIO  -  0.9350 
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MODEL  =  BACKFIRE 
ALT  =  !0  KM 
HEADING  =  0.0 

YAW  =  PITCH  =  ROLL  =  0.0  DEG 
EARTH  CENTRAL  ANG  =  0.0  DEG 
TGT  TEMP*  250  K 
SENSOR  ALT  =  0  K 
CENTRAL  OBS  FREQ  =  35  GHZ 
BANDWIDTH  =  1  GHZ 
ATMO  *  US  STD 
SOLAR  IRRAD  -  0.0  W  /  M**2 


OBSERVER  DISTANCE  *  10  KM 
MEAN  ATMO  XM  IS  SION  *  0.88 
TOT  PROJ  AREA  =  290.9  SQ  M 
BACKGND  TEMP  =  36.3  K 
MEAN  TGT  TEMP  =  252.7  K 
MEAN  APPARENT  TEMP  =  190.5  K 
ME.  M  INHERENT  TEMP  =  216.4 
TGT  APPARENT  CONTRAST  =  5.244 
TGT  INHERENT  CONTRAST  *  5.959 
APPARENT  /  INHERENT  RATIO  =  0.88 


SURFACE  =  LAND 
SFC  TEMP  =  290  K 
SFC  TYPE  =  VEGETATION 
SFC  FRACTIONAL  COVER  *  0  75 
SFC  DIELECTRIC  CONST  =  10 
SFC  BACKS CATT  COEF  =  0.2 
CLOUD  =  WATER 
CLOUD  BOTTOM  =  1  KM 
CLOUD  THICK  =  I  KM 
CLOUD  LWC  -  0.25  GM  /  M**3 


MODEL  =  BACKFIRE 
ALT  =  10  KM 
HEADING  =  0.0 

YAW  =  PITCH  =  ROLL  =0.0  DEG 
EARTH  CENTRAL  ANG  =  0.0  DEG 
TGT  TEMP  =  250  K 
SENSOR  ALT  =  20  KM 
CENTRAL  OBS  FREQ  =  35  GHZ 
BANDWIDTH  =  1  GHZ 
ATMO  »  US  STD 
SOLAR  IRRAD  «  0.0  W  M**2 


SURFACE  =  LAND 
SFCTEMP  =  290  K 
SFC  TYPE  =  VEGETATION 
SFC  FRACT  IONAL  COVER  =0.75 
SFC  DIELECTRIC  CONST  =  10 
SFC  BACKSCATT  COEF  =  0.2 
CLOUD  =  NONE 
CLOUD  BOTTOM  = 

CLOUD  THICK  = 

CLOUD  LWC  = 


OBSERVER  DISTANCE  =  10  KM 
MEAN  ATMO  XM1SSION  =  0.9963 
TGT  PROJ  AREA  =  288.3  SQ  M 
BACKGND  TEMP  =  260.9  K 
MEAN  TGT  TEMP  =  30.96  K 
MEAN  APPARENT  TEMP  =  229.1  K 
MEAN  INHERENT  TEMP  =  229.9  K 
TGT  APPARENT  CONTRAST  =  -0.8781 
TGT  INHERENT  CONTRAST  =  -0.88 1 3 
APPARENT  /  INHERENT  RATIO  =  0.9963 
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THE  SUN  BECOMES  A  MORE  IMPORTANT  GLINT  SOURCE 
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MODELS  FOR  WARM  CLOUD  TYPES,  ICE  CLOUDS,  AND  PRECIPITATION 
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ABSORPTION  CROSS  SECTION  MEASUREMENTS  OF  SCHUMANN-RUNGE 
CONTINUUM  OF  02  AT  78  K  AND  295  K 

K.  Yoshino,  J.R.  Esmond,  W.H.  Parkinson 
Harvard-Smithsonian  Center  for  Astrophysics,  Cambridge,  MA  02138 

We  measured  the  photoabsorpdon  cross  sections  of  the  S-R  continuum  at  78  K  and  295  K 
in  the  wavelength  region  130-175  nm.  The  1-m  Seya-Namioka  vacuum  spectrometer  on  the  BL- 
12A  beam  line  at  die  Photon  Factory,  KEK,  Japan  was  used,  in  die  first  order,  with  synchrotron 
radiation  as  a  background  continuum.  Die  measurements  have  been  performed  on  a  relative  scale 
in  scanning  mode  and  also  on  the  absolute  scale  at  every  5  nm  interval.  These  absolute  cross 
sections  measurements  of  O2  have  been  used  to  put  the  relative  cross  sections  on  a  firm  absolute 
basis  throughout  die  region  130-175  nm.  Temperature  effects  on  the  cross  sections  of  die  S-R 
continuum  are  clearly  seen  and  cross  sections  at  78  K  are  about  80%  of  those  at  295  K.  The 
observed  cross  sections  at  78  K  are  limited  to  absorption  originating  from  only  the  v"=0  level. 

This  work  is  supported  by  NASA  grant  NAG5-484  to  Smithsonian  Astrophysical 
Observatory. 


Absorption  Cross  Section  Measurements  of  Schumann-Runge 
Continuum  of  Oj  at  78  K  and  295  K 
K.  Yoshino,  J.R.  Esmond  and  W.H.  Parkinson 
Harvard-Smithsonian  Center  for  Astrophysics,  Cambridge,  MA  0S1S8,  U.S.A. 


We  measured  the  photoabsorption  cross  sections  of  the  S-R  continuum  at  78  K  and  295 
K  in  the  wavelength  region  130-175  nm.  The  1-m  Seya-Namioka  vacuum  spectrometer  on 
the  BL-12A  beam  line  at  the  Photon  Factory,  KEK,  Japan  was  used,  in  the  first  order,  with 
synchrotron  radiation  as  a  background  continuum.  The  measurements  have  been  performed 
on  a  relative  scale  in  scanning  mode  and  also  on  the  absolute  scale  at  every  5  nm  interval. 
These  absolute  cross  sections  measurements  of  02  have  been  used  to  put  the  relative  cross 
sections  on  a  firm  absolute  basis  throughout  the  region  130-175  nm.  Temperature  effects  on 
the  cross  sections  of  the  S-R  continuum  are  clearly  seen  and  cross  sections  at  78  K  are 
about  80%  of  those  at  295  K.  The  observed  cross  sections  at  78  K  are  limited  to  absorption 
originating  from  only  the  v'=0  level. 

This  work  is  supported  by  NASA  grant  NAG5-484  to  Smithsonian  Astrophysical 
Observatory. 
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Absolute  Absorption  Cross  Sections  Measurements 


59.8  239. 


section  measurements  of  Schumann-Runge 
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THE  SPECTRAL  RADIATION  EXPERIMENT  (SPECTRE):  AN  OVERVIEW  - 
CLEAR-SKY  OBSERVATIONS  AND  THEIR  USE  IN  ICRCCM 

R.G.  Ellingson 

Department  of  Meteorology,  University  of  Maryland,  College  Park,  MD  20742 

WJ.  Wiscombe,  V.  Kunde,  H.  Melfi 
NASA  Goddard  Space  Flight  Center,  Code  613,  Gteenbelt,  MD  20771 

J.  DeLuisi 

Air  Resources  Laboratory,  NOAA/Environmental  Research  Laboratory,  Boulder,  CO  80303 

D.Murcray 

Department  of  Physics,  University  of  Denver,  2400  S.  Vine,  Denver,  CO  80208 

W.  Smith 

Space  Science  &  Engineering  Center,  University  of  Wisconsin, 

1225  W.  Dayton  Street,  Madison,  WI  53706 

Surface-based,  high  resolution  observations  of  die  spectral  radiance  in  die  3  to  20  pm  region 
were  obtained  simultaneously  by  three  F11K  spectrometers  concurrently  with  insitu  and  remote 
observations  of  atmospheric  temperature  and  water  vapor  profiles  and  trace  gases  CO2,  CH4, 

N2O,  O3,  FI  1  and  F12  as  pan  of  the  FIRE  Gnus  II  field  experiment  The  spectrometers  were 
self-  and  inter-calibrated  in  the  field  to  achieve  close  to  1%  accuracy.  The  presentation  will 
summarize  the  data  obtained  for  selected  clear-sky  cases,  and  comparisons  of  model  calculations 
with  observations  will  be  used  to  illustrate  how  the  data  may  be  used  to  clarify  uncertainties 
associated  with  calculations  of  the  effects  of  trace  gases  as  round  by  die  international 
Intercomparison  of  Radiation  Codes  in  Climate  Models  (ICRCCM). 
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ON  THE  USE  OF  GROUND-BASED  REMOTE  SENSORS  TO  PROVIDE 
TEMPERATURE  AND  MOISTURE  SOUNDINGS  FOR  ATMOSPHERIC 

TRANSMISSION  MODELS 


J.C.  Liljegren 

Battcllc,  Pacific  Northwest  Laboratory,  Richland,  WA  99352 

R.O.  Knuteson 

University  of Wisconsin,  1225  W.  Dayton  Street,  Madison,  WI,  53706 


A  general  goal  of  die  Atmospheric  Radiation  Measurement  (ARM)  Program  of  the  DOE  is  to 
advance  the  state-of-the-art  for  models  of  atmospheric  transmission  through  comparison  of  model- 
calculated  radiance  spectra  with  measured  spectra.  ARM  will  rely  primarily  on  ground-based 
remote  sensors  to  acquire  the  input  temperature  and  moisture  profiles  simultaneously  with  the 
radiance  measurements.  To  support  this  strategy,  we  present  comparisons  of  measured  infrared 
radiance  spectra  with  spectra  calculated  using  FASCQDE  for  two  balloon-borne  systems,  NCAR- 
GLASS  and  NWS,  and  two  remote  sensing  systems,  die  NOAA/WPL  6-channel  microwave 
profiler  (MWP)  and  the  High-resolution  Infrared  Sounder  (HIS),  which  were  co-located  and 
simultaneously  operated  dining  the  Ground-based  Atmospheric  Profiling  Experiment  (GAPEX) 
described  by  Smith  etaL  (1990).  Our  results  show  that,  at  least  for  die  relatively  dry  conditions 
which  prevailed  dnring  GAPEX,  the  calculations  based  on  the  remote  sensor-based  profiles  were 
at  least  as  accurate  as  those  using  the  balloon-borne  systems. 


On  the  Use  of  Ground-Based  Remote  Sensing 
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Radiance  (HIS) 


The  GAPEX  Data  Set 

i-based  Atmospheric  Profiling  Experiment 

Denver  Stapleton  Airport,  29  Oct  -  4  Nov.  1988 
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Retrieved  vertical  profiles  of  temperature  and  humidity 


Testing  Method 
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GAPEX  31  Oct  88 
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MODEL  CLIMATOLOGIES  OF  TRACE  SPECIES  IN  THE  ATMOSPHERE 

M.E.  Summers 

Naval  Research  Laboratory,  Upper  Atmospheric  Physics  Branch,  Washington,  DC  20375 

WJ.  Sawchuck 

Computational  Physics,  Inc.,  P.O.  Box  788,  Annandale,  VA  22003 

GJP.  Anderson 

Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  AFB,  MA  01731 

A  set  of  model  climatologies  of  important  chemically  and  optically  active  trace  species  in  the 
atmosphere  has  been  developed  for  use  with  atmospheric  radianceAransmittance  codes.  The  new 
climatologies  indude  monthly  mean  pressure/temperature/density  fields  along  with  N2, 02, 03, 0, 
H2O,  GH4,  CO2,  CO,  NQ2.  N2O,  and  HNO3  winch  represent  zonally  averaged  abundances  for 
the  altitude  range  0  - 120  lam,  at  10  degree  latitude  intervals.  Fra  O3  and  O,  day/night  ratios  are 
tabulated.  The  climatologies  incorporate  observations  of  trace  species  abundances  where  available 
and  validated  modd  extrapolations  in  other  regions.  Variances  about  the  monthly  mean  values 
have  been  deduced  from  selected  observational  data  sets  and  provide  estimates  of  variability  in  the 
constituent  fields.  The  generation  of  die  climatologies  will  be  discussed  and  selected  portions  will 
be  shown.  Community  access  to  the  climatologies  will  also  be  discussed. 


MODEL  CLIMATOLOGIES  OF  TRACE  SPECIES  IN  THE  ATMOSPHERE 


M.E.  Suntn,  Naval  Research  Laboratory 
W.J.  Savchuck,  Computational  Physics,  Inc. 

G.P.  Anderson,  Geophysics  Directorate,  Phillips  Laboratory 


Presentation  at: 

ANNUAL  REVIEW  CONFERENCE  ON  ATMOSPHERIC  TRANSMISSION  MODELS 
2-3  June  1992 

Geophysics  Directorate,  Phillips  Laboratory 
Hanscoa  AFB,  MA  01742 
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MODEL  CLIMATOLOGIES  OF  TRACE  SPECIES  IN  THE  ATMOSPHERE 


M.E.  Summers 

Naval  Research  Laboratory,  Upper  Atmospheric  Physics  Branch,  Washington,  DC  20375 

WJ.  Sawchuck 

Computational  Physics,  Inc.,  P.O.  Box  788,  Annandale,  VA  22003 
G.P.  Anderson 

Geophysics  Directorate,  Phillips  Laboratory,  Hanscom  AFB,  MA  01731 
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A  set  of  model  climatologies  of  important  chemically  and  optically  active  trace 
species  in  the  atmosphere  has  been  developed  for  use  with  atmospheric 
radiancc/transmittance  codes.  The  new  climatologies  include  monthly  mean 
pressure/temperature/density  fields  along  with  N2, 02, 03. 0,  H2O,  CH4,  CO2,  CO,  NO2, 
N2O,  and  HNO3  which  represent  zonally  averaged  abundances  for  the  altitude  range  0  - 
120  km,  at  10  degree  latitude  intervals.  For  O3  and  O,  day/night  ratios  are  tabulated.  The 
climatologies  incorporate  observations  of  trace  species  abundances  where  available  and 
validated  model  extrapolations  in  other  regions.  Variances  about  the  monthly  mean  values 
have  been  deduced  from  selected  observational  data  sets  and  provide  estimates  of  variability 
in  the  constituent  fields.  The  generation  of  the  climatologies  will  be  discussed  and  selected 
portions  will  be  shown.  Community  access  to  the  climatologies  will  also  be  discussed. 
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I.  INTRODUCTION 
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to  generate  this  climatological  database  is  provided  in  a  separate  report. 


r  i 


L  A 


II.  DATABASE  CONTENT 


The  Trace  Constituent  Climatology  Database  currently  consists  of  27  individual  pa¬ 
rameter  fields.  These  fields  are  classified  as  follows: 

Reference  Atmosphere  Climatology 

(1)  Pressure 

(2)  Temperature 

(3)  Density 

Ozone  Diurnal  Variability 

(4)  Average  Night/Day  Variation  of  O3  mixing  ratios 

Trace  Species  Climatology  (Mixing  Ratios) 

(5)  O, 

(6)  H20 

(7)  CO 

(8)  CO, 

(9)  CHA 

(10)  N02 

(11)  n2o 

(12)  HNO3 

(13)  O 

(14)  O, 

(15)  N2 

Trace  Species  Variability  (Std.  Dev.) 

(16-17)  03 
(18-19)  H%0 
(20-21)  NOi 
(22-23)  HN03 
(24-25)  N20 
(26-27)  CHa 
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25  km  -  120  km  (5  km  resolution) 


mbar):  80  km 


tr  AmniAwnronwv  dec 
TIKE  (month) 


MIXING  RATIO  (ppbv):  JAN  N,0  MIXING  RATIO 


UTtTUDt 


bv):  JAN  _ HNO.  MIXING  RATIO  (ppbv):  FEB 


NICHT/DAY  RATIO:  JAN  0,  NICHT/DAY  RATIO 
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FAR- WING  LINE  SHAPE  CONTRIBUTION  TO  THE  WATER  CONTINUUM: 

RECENT  RESULTS  AND  FUTURE  DIRECTIONS 

R.H.  Tipping 

Department  of  Physics  and  Astronomy,  University  of  Alabama,  Tuscaloosa,  AL  35487 

Q.Ma 

Institute  for  Space  Studies,  Goddard  Space  Flight  Center,  New  York,  NY  10025 

A  far-wing  line  shape  theory  has  been  developed  and  applied  to  die  calculation  of  both  self- 
and  N^broadened  water  lines.  This  theory  is  based  on  the  quasistatic  and  binary  collision 
approximation  and  utilizes  only  known  parameters  obtainable  from  other  measurements  (line 
strengths,  positions,  interaction  potential  parameters,  and  multipole  moments).  The  resulting 

•  continuous  absorption  coefficient  is  in  good  agreement,  in  both  its  magnitude  and  its  temperature 
dependence,  with  existing  experimental  laboratory  measurements.  Some  refinements  in  the  model 
ana  future  directions  will  be  discussed  briefly. 
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INTRODUCTION 

1.1)  IMPORTANCE  OF  MILLIMETER  AND 

INFRARED  RADIATION  IN  ATMOSPHERES 

(1)  Controls  the  energy  balance  and 
affects  the  climate  of  the  Earth. 

(2)  Widely  used  in  remote  sensing, 
atmospheric  communications,  etc. 

(3)  To  determine  characteristics  of 
planetary  atmospheres  from 
experiment  data. 


r  ^ 


1.2)  ATMOSPHERIC  WINDOWS  AND 
CONTINUUM  ABSORPTION 

OPAQUE  REGIONS  AND  WINDOWS: 

The  opaque  regions  are  mainly  determined 
by  H20  lines,  along  with  02  and  C02 
lines.  Between  them  there  are  ”WlNDOWSH. 

THE  WEAK  ABSORPTION  IN  WINDOWS: 

(1)  Some  weak  local  lines  of  H20,  02 
and  C02« 

(2)  A  gradually  varying  background, 
i.e.  continuum  absorption. 
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1.3)  BASIC  CHARACTERISTICS  OF  CONTINUA 

There  are  two  kinds  of  water  continua: 

(1)  H2o  -  H20  (self)  continuum 

(2)  H20  -  N2,  02  (foreign)  continuum 
Usually,  the  self -continuum  is  dominant 

CHACTERISTICS  OF  THE  SELF-CONTINUUM : 

(1)  Varies  smoothly  with  frequency. 

(2)  Proportional  to  n2. 

(3)  Very  strong  negative  temperature 
dependence • 


1.4)  MAIN  THEORIES  OF  WATER  CONTINUUM 

THREE  MECHANISMS  FOR  SELF-CONTINUUM: 

(1)  Collision  broadened  far-wings  of 
allowed  water  vapor  lines. 

(2)  Collision-induced  absorption. 

( 3 )  Water  dimers • 

All  three  mechanisms  vary  as  n2, 
but  there  are  large  differences  in  the 
temperature  dependencies  and  in  the 
magnitudes  of  the  absorption. 
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2.  FAR-WING  LIME  SHAPE  THEORY 

a(w)  =  (4trJ/3^c)n,u  tanh(flu/2kT)  [F(u)  +  F(-w)  ]  (1) 

00 

F (u)  =  Im  (2trJi)‘‘  $  du'  Tr[M»*  (u-u'-L,)'1 

-00 

X  <(u'-Lt)1  p«“»>k  p«M.]  (2) 
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APPROXIMATIONS 

1.  BINARY  COLLISION  APPROXIMATION 

p<"»  =  n„  e-v<"/ta 
v(r)  =  V„110(r)  +  V,„(r) 

V^Cr)  =  Co'6 [  ( ct / r ) 1  +  (<J/r)6] 

2 •  QUASISTATIC.  APPROXIMATION 
<(w'-Ll)'1  p(inl)>b 
<...>->  47t  J  ...  r2dr 
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3.  BAND-INDEPENDENT  APPROXIMATION 


a(oj)  =  E  avivj(w) 

vivj 


ayjv)(cj)  =  n,  E  Sv1¥j[<j  sinh(/jw/2kT)  ]/ [cjvivj  sinh(^wvivj/2kT)  ] 


Uv.vjX) 


x  ^  /  (cj-ajvivj) 2  +  X.(cj+wvivj)  /  (w+wvWj)2> 


Y  ^ 
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INPUT  DATA 

POTENTIAL  PARAMETERS 
C,  a,  t,  H,  Q 

SPECTROSCOPIC  PARAMETERS 


(HITRAN  DATABASE) 


THE  LINE  SHAPE  FUNCTIONS,  /,(<*<)  *nd  /-(w)  (c®'1/**®) 


|.QE*QI 


I.QE+CO 


1  - QE-C !  == 


1.0E-G2  = 


1  .UE-03=: 
tz 


1.0E-CM 


10000 


FIG.  (1) 


FREQUENCY  « (cm •') 
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ABSORPTION  COEFFICIENT  a(u)  (cm1  molecule'1  elm*') 


f  ^ 
L  A 


f  ^ 

L  A 


"h  Suaci-'  a-t 

0  —  1 
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FREQUENCY  w  (cm'1) 


ABSORPTION  COEFFICIENT  o(w)  (cm1  molecule*'  atm* ) 


<T.(w)  (cm'V»tm) 


ABSORPTION  COEFFICIENT  n(u>)  (cm*  molecule*1  »tm*') 


:ule*'  »tm* 


UtO  (  -1  V_VA  V*.  V-k 


WAVFNUMRFR 


180 
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4.  FUTURE  DIRECTIONS 

A.  MORE  REALISTIC  POTENTIAL  MODEL 

B.  CORRECTIONS  TO  QUASISTATIC  APPROXIMATION 

C.  INTERFERENCE  BETWEEN  ALLOWED  AND 
INDUCED  DIPOLE  MOMENTS 

D.  VALIDATION  WITH  NEW  DATA  WHEN  AVAILABLE 
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TOWARD  AN  IMPROVED  WATER  VAPOR  CONTINUUM  MODEL 

S.A.  Clough 

Atmospheric  and  Environmental  Research,  Inc.,  840  Manorial  Drive,  Cambridge,  MA  02139 

A  model  for  die  water  vapor  continuum  based  on  the  measurements  of  Burch  [1],  has  been 
reported  by  Clough  et  al.  [2,3].  Since  the  definition  of  the  continuum  is  inextricably  related  to  the 
line  contribution  associated  with  die  radiative  transfer  program  with  which  it  is  used,  a  consistent 
analysis  of  the  data  is  essential  to  die  proper  treatment  of  die  problem.  To  obtain  an  improved 
water  vapor  continuum  model,  the  load  line  contribution  consistent  with  the  line-by-line  program 
FASGODE  has  been  carefully  and  systematically  removed  from  the  extinction  measurements  made 
in  the  microwindows  of  the  water  vapor  spectrum  by  Burch  [1].  Other  measurements  based  on  the 
transmittance  over  a  homogeneous  atmospheric  path  and  radiances  for  an  inhomogeneous 
atmospheric  path  are  being  considered  in  the  analysis.  An  improved  line  shape  including  the 
effects  of  detail  balance  is  being  developed  on  the  basis  of  these  measurements.  The  results  from 
the  present  model  are  compared  with  recent  theoretical  results  of  Ma  and  Tipping  [4]. 

1.  Bwch.  D-  E„  Continuum  absorption  by  H?Q.  Rep.  AFGL-TR-8 1 -0300. 32  pp.  Air  Force 
Geopbys.  Lab.,  Hanscom  AFB,  MA,  1981.  ADA1 12264 

2.  Clough,  S.  A.,  F.  X.  Kneizys,  R.  W.  Davies,  R.  Gamache,  and  R.  H.  Tipping,  Theoretical 
line  shape  for  H2O  vapor  application  to  the  continuum-  in  Atmospheric  Water  Vapor,  edited  bv 
A.  Deepak,  T.  D.Wilkerson,  and  L.  H.  Rhunke,  pp.  25-46,  Academic  Press,  London,  1980. 

3.  Clough,  S.  A.,  F.  X.  Kneizys,  and  R.  W.  Davies,  Line  shape  and  the  water  vapor 
continuum.  Atmospheric  Research,  22, 229-241, 1989b. 

4.  Ma,  Q.,  and  R.  H.  Tipping,  A  far  wing  line  shape  theory  and  its  application  to  the  water 
continuum  absorption  in  die  infrared  region  (II),  J.  Chem.  Phvs..  in  press,  1992. 
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A  CORRELATED  K-DISTRIBUTION  MODEL  OF  THE  HEATING  RATES  FOR 
CH4  AND  N20  IN  THE  ATMOSPHERE  BETWEEN  0  AND  60KM* 


K.  Grant,  A.S.  Gross 

Lawrence  Livermore  National  Laboratory,  P.O.  Box  808,  L-262,  Livermore,  CA  94551-9900 


An  infrared  radiative  transfer  model  using  die  correlated-k  technique  is  being  developed  for 
use  in  atmospheric  chemistry  and  climate  models  at  LLNL,  as  well  as  for  stand-alone  radiative 
diagnostics.  This  implementation  of  die  correlated-k  algorithm  will  provide  a  computationally 
efficient  means  to  calculate  atmospheric  adiative  fluxes  and  heating  rates.  It  will  also  allow  for 
both  time-varying  gas  concentrations  and  far  the  rapid  inclusion  of  new  gases  into  the  model 
We  present  the  results  of  our  prototype  model  of  absorption  by  CH4  and  N2O  in  the  7-9 
micron  wavelength  band.  The  correlated-k  algorithm  maps  the  highly  structured  absorption 
coefficient  versus  wavelength  relationship  into  a  substantially  smoother  monotonic  relationship 
between  die  absorption  coefficient  and  its  cumulative  probability  distribution.  This  provides  a 
substantial  reduction  in  computational  complexity  relative  to  line-by-line  calculations,  while 
maintaining  a  high  degree  of  accuracy. 


♦Work  performed  under  the  auspices  of  the  U.S.  Dept  of  Energy  by  the  Lawrence  Livermore 
National  Laboratory  under  contract  No.  W-7405-Eng-48. 
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Abstract 

An  infrared  radiative  transfer  model  using  the  correlated-*  technique  is  being  developed  for  use 
in  atmospheric  chemistry  and  climate  models  at  LLNL,  as  well  as  for  stand-alone  radiative 
diagnostics.  This  implementation  of  the  correlated-*  algorithm  will  provide  a  computationally 
efficient  means  to  calculate  atmospheric  radiative  fluxes  and  beating  rates.  It  will  also  allow  for 
both  time-varying  gas  concentrations  and  for  the  rapid  inclusion  of  new  gases  into  the  model 

We  present  both  a  general  background  of  the  c-k  method  and  the  results  of  applying  our 
prototype  model  to  absorption  by  CH4  and  NjO  in  the  7-9  micron  wavelength  band.  The 
correlated-*  algorithm  maps  the  highly  structured  absorption  coefficient  versus  wavelength 
relationship  into  a  substantially  smoother  monotonic  relationship  between  the  absorption 
coefficient  and  its  cumulative  probability  distribution.  This  provides  a  substantial  reduction  in 
computational  complexity,  relative  to  line-by-line  calculations,  while  maintaining  a  high  degree 
of  accuracy. 

L  Introduction 

Infrared  radiative  transfer  models  are  used  at  LLNL  both  as  stand-alone  diagnostic  tools 
and  as  important  components  of  more  general  atmospheric  models.  As  diagnostic  models,  we  use 
radiative  transfer  models  to  calculate  the  radiative  forcing  to  the  surface-troposphere  system  due 
to  potential  changes  in  trace  gases  concentrations,  aerosols,  or  clouds.  For  example,  such 
radiative  forcing  calculations  are  an  essential  step  in  deriving  global  wanning  potentials  (GWPs) 
for  trace  gases.  These  GWPs  provide  an  approximate  index  of  the  time-cumulative  radiative 
effects  of  a  unit  emission  of  a  specified  trace  gas  relative  to  the  comparable  effect  for  CO? 
Another  diagnostic  application  of  our  efforts  in  atmospheric  radiative  transfer  modeling  has  been 
the  calculation  of  stratospheric  temperature  changes  forced  by  observed  changes  of 

ozone  mixing  ratio  profiles. 

As  integral  parts  of  global  and  regional  chemistry  mod  .  md  of  climate  models,  we  use 
infrared  radiative  transfer  models  in  calculating  vertical  profiles  of  net  heating  rates.  Such 
heating  rate  profiles  are  a  vital  part  both  of  calculating  the  atmospheric  circulation  within  global 
models  and  of  estimating  the  chemical  and  temperature  feedback  effects  from  potential 
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anthropogenic  trace  gas  emissions.  As  a  requirement  of  our  global  modeling  efforts,  we  are 
emphasizing  the  development  of  radiative  transfer  models  that  maintain  accuracy  into  the  upper 
stratosphere  arid  that  are  implemented  with  the  flexibility  to  handle  time  varying  profiles  of 
temperature  and  radiatively  active  overlapping  trace  gases. 

A  current  focus  of  our  efforts  is  the  development  of  an  infrared  radiative  transfer  model 
which  can  provide  a  substantial  reduction  in  computational  complexity,  relative  to  line  by  line 
calculations,  while  maintaining  a  high  degree  of  accuracy.  In  a  recent  publication  (Ladas  and 
Oinas,  1991,  “LO"),  a  method  was  described  for  treating  transmission  calculations  in  infrared 
radiative  transfer  problems  for  a  vertically  inhomogeneous  atmosphere  in  which  there  is  non¬ 
gray  gaseous  absorption.  The  method  known  as  the  correlated  £ -distribution  utilizes  a  mapping 
of  die  opacity-frequency  relation  into  an  opacity-probability  relation  within  a  particular 
frequency  interval.  The  probability  variable  g(k),  the  cumulative  distribution  function,  is  defined 


g(k) =]/&)<*  • 


where  flk')dk'  is  the  fraction  of  the  frequency  interval  occupied  by  absorption  coefficients 
between  k'  and  k'+dk‘  (Goody  and  Jung,  1989,  lX31";  Goody  a  al.,  1989,  “G2”;  and  West  er  al., 
1990,  “Wl").  The  limits  of  g(k)  range  between  0  and  1  within  the  frequency  interval  The  inverse 
of  Equation  1,  k(g),  the  ^-distribution,  has  been  shown  by  LO,  G2,  and  Wl  to  be  a  monotonic 
function  across  the  frequency  interval  for  a  particular  atmospheric  layer.  The  correlated-k 
method  can  mathematically  provide  an  exact  procedure  for  calculating  the  transmission  and 
heating  rates  in  a  homogeneous  atmosphere.  For  the  case  of  inhomogeneous  atmospheric  paths 
the  method  is,  in  practice,  inherently  inexact  since  somewhat  different  sets  of  frequencies  will 
associate  with  a  given  ordering  of  the  k  terms  as  the  pressure  and  temperature  vary  over  the  path. 
If  the  fractional  error  caused  by  the  use  of  the  correlated-k  method  to  calculate  the  transmission 
for  inhomogeneous  paths  is  tolerable,  then  the  method  offers  an  attractive  means  for  computing 
radiation  fields  in  real  atmospheres. 

The  calculation  of  the  transmission  can  be  expressed  in  the  three  physically  equivalent 
forms: 


T(u)  =  11  dvf  exp  (~k^u)dv  , 

-  J/WexpftfiOd?  , 

o 

t 

=  Jexp(-*(g)u)dg  , 


where  u  is  the  absorber  column  density.  Using  the  £ -distribution  form,  the  calculation  can  be 
performed  using  far  fewer  k-g  points  than  the  same  calculation  using  k-v  (frequency)  points. 
Thus  the  £ -distribution  method  has  the  potential  of  being  a  much  more  computationally  efficient 
method  of  doing  radiative  transfer  calculations  within  the  earth’s  atmosphere  when  compared  to 
line  by  line  or  narrow  band  models  while  retaining  many  of  the  numerical  advantages  of  an 
explicitly  exponential,  monochromatic  calculation.  For  example,  scattering  may  be  readily 
handled  with  ti».  correlated  k -distribution  model,  but  not  with  band  models. 
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The  opacity  model  used  by  LO  for  the  molecules  CH*  and  N2O  was  a  Malkmus  band 
model  (Malkmus,  1967;  Rodgers,  1968)  for  the  transmission  function, 

T(u)  =  exp|-  *  fl/2[(l  +  4Su/xB)m  - 1]|  ,  ^ 

where  the  parameters  S  and  B  are  obtained  by  fits  to  the  mean  absorption  line  strength,  mean 
absorption  line  half  width,  and  mean  line  spacing  in  the  frequency  interval.  The  probability 
function /ft)  in  Equation  1  has  been  shown  to  be  the  inverse  Laplace  transform  of  Equation  3  by 
Domoto  (1974)  and  is  given  by  the  expression  (LO); 

/(kJsO-SiT^CSS^expIttB/d^-Syk-t/S]} .  (4) 

The  ^-distributions  used  by  LO  for  CH*  and  N2O  were  obtained  using  the/?*)  relations  derived 
from  the  Malkmus  band  model.  Excellent  fits  to  the  CH4  and  N2O  heating  rates  were  obtained  by 
LO  using  the  Malkmus  model  with  parameters  empirically  fit  to  give  the  best  agreement  with 
line  by  line  calculations. 

An  more  direct  and  inherently  more  accurate  method  of  obtaining  *-distributions  has  been 
used  by  G2  and  by  Freedman  and  Pollack  (1990,  “FP”).  This  method  involves  a  direct  line  by 
line  calculation  of  the  opacity-frequency  relation  and  a  direct  numerical  sorting  of  this  data  to 
obtain/?*)  and  the  *-distribution,  *?g).  The  advantage  of  doing  this  is  that  the  Malkmus 
formulation  is  bypassed  and  the  assumptions,  such  as  Loren tz  line  stupe  introduced  into  the 
problem  by  the  Malkmus  model  are  avoided. 

The  main  purpose  of  this  paper  is  to  review  the  theory  necessary  to  implement  the 
correlated-*  model  for  a  vertically  inhomogeneous  atmosphere,  and  to  apply  the  direct  FP 
method  to  calculate  *-distributions  for  the  transmission  of  the  molecules  CH4  and  N2O.  For  the 
wavelength  interval  of  7  to  9  microns  and  altitude  range  of  0  to  62  km,  a  calculation  of  the 
atmospheric  cooling  rates,  upward  fluxes,  and  downward  fluxes  will  be  made  for  each  molecule 
and  compared  to  the  LO  values. 

n.  The  K-Distribution  Method 


For  a  non-homogeneous  atmospheric  path,  containing  j  discrete  layers,  the  transmission 
expression  given  by  Equation  2  becomes; 


T[ii )  =  1/dv 


av  V.  /  / 

=  J'«P (“X**/4")  ■ 

0  V  /  / 


dv  , 


(5a) 

(5b) 


The  exact  equivalence,  of  Equations  Sa  and  3b  for  an  inhomogeneous  atmosphere  is  limited 
to  a  few  special  cases.  G1  has  shown  that  these  expressions  are  equivalent  for  the  case  of 
*«*  =*«©/)  .where  the  subscript  “o” ,  indicates  a  particular  layer  and  /.  is  a  function  of  the 
pressure  and  temperature  of  the  /  *  layer.  G2  has  shown  that  the  two  expressions  are  also 
equivalent  for  die  case  of  the  strong  line  limit  for  pressure  broadened  lines.  When  the  exponential 
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terms  in  Equations  Sa  and  Sb  are  expanded  and  the  integrals  are  compared  term  by  term,  GI  has 
shown  that  in  the  weak  line  limit. 


the  two  expressions  are  equivalent  Gl  also  states  that  since  terms  of  the  kind. 


contain  the  expressions  (k.AUj)  which  are  equivalent  in  the  two  transmission  models,  the 
k  distribution  method  will  oc  a  more  accurate  method  than  the  weak  line  approximation.  In 
general,  the  term  integrals  which  involve  the  cross  product  terms  will  not  be  equivalent  unless 
there  is  a  unique  k(g)  relation  which  is  correlated  for  all  the  layers.  Both  LO  and  W1  have 
discussed  the  validity  of  the  /^-distribution  method  from  this  point  of  view.  For  the  method  to 
work  precisely,  the  k(g)  relations  must  be  such  that  a  particular  g  value  selects  the  same  set  of 
frequencies  in  all  the  layers;  i.e.  the  frequency  ordering  of  the  opacity  groups  is  preserved  layer 
to  layer.  In  LO  and  W1  examples  are  presented  which  show  the  consequences  that  can  occur 
when  specific  g  values  do  not  select  the  same  frequency  sets  in  different  layers.  Via  the  separate 
k(g)  relationships  for  two  layers,  the  k(  v)  distribution  for  a  layer  can  be  scaled  on  to  the  range  of 
k  values  appropriate  for  the  second  layer  by  aeleeting  a  set  of  frequencies  appropriate  to  a 
particular  g  value  for  the  first  layer  and  using  this  set  of  frequencies  with  the  k  value,  for  the 
same  g,  in  the  second  layer  as  the  scaling  relation.  The  scaled  k(v)  relation,  thus  obtained  for  the 
second  layer,  will  differ  from  the  original  k(v)  relation  for  the  second  layer,  with  the  difference 
increasing  as  the  correlation  between  the  k{v)  relations  of  the  two  layers  decreases.  The  k- 
distribution  derived  from  the  scaled  k(v)  relation  will  become  significantly  non  monotonic  as 
the  correlation  between  the  layers  decreases;  Le.  the  k  value  for  a  given  g  will  no  longer  show  a 
small  dispersion  about  a  well  defined  mean  value.  As  pointed  out  in  G2,  the  ^-distribution 
method  cannot  be  valid  for  all  conditions  and  its  use  under  all  circumstances  mna  be 
accompanied  by  numerical  validity  tests  or  comparison  with  asymptotic  forms.  The  monotonic 
nature  of  the  k-distributions  allow  transmission  calculations  to  be  performed,  over  a  frequency 
interval,  using  a  few  tens  of  k(g)  points  as  opposed  to  large  numbers  of  k(v)  points  (Ao-Auh/4; 
AUH-line  half  width)  required  for  a  line  by  line  calculation.  The  resulting  computational 
economy  using  die  k-distribution  method  is  sufficient  to  warrant  the  numerical  validity  tests 
necessary  for  its  use. 

The  calculation  of  the  k-distributions  for  CH«  and  NjO  by  the  FP  method  contains  the 
following  steps.  First,  the  HTTRAN-91  database  (Rothman  a  al„  1991)  is  utilized  to  determine 
the  line  transitions  and  physical  properties  of  the  selected  lines.  Second,  a  version  of  the 
FASCODE2  code  (Qough  «  at.,  1986),  modified  by  FP,  is  used  to  calculate  a  finely  gridded 
(Ao~Auh^)  aet  of  monochromatic  absorption  coefficients,  with  full  allowance  for  the  overlap  of 
neighboring  lines,  for  each  layer  in  the  atmosphere.  Third,  a  sorting  code  developed  by  FP, 
ASSORT,  is  used  to  calculate  th cflk),  g(k),  and  k(g)  functions  for  each  homogeneous  layer.  A 
block  diagram  of  the  process  is  shown  in  Figure  1.  Initial  input  to  the  SELECT  code  which 
accesses  the  HTTRAN-9 1  database  consists  of  the  wave  number  range,  and  the  choice  of 
molecular  species  and  isotopes.  The  SELECT  code  output  contains  the  central  wave  number, 
intensity,  half  width,  and  lower  state  energy  for  each  line  as  well  as  an  intensity  weighted 
average  half  width  for  all  the  lines.  The  weighted  average  half  width  calculation  in  the  SELECT 
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code  wis  a  modification  added  by  FP.  The  modified  FASCODE  program  takes  the  line  data 
generated  by  the  SELECT  code,  fits  an  absorption  line  profile  (Voigt  Profile)  to  each  line  and 
calculates  the  absorption  coefficient  k(cm2/air  mol)  as  a  function  of  wave  number.  The  wave 
number  grid  spacing  is  controlled  by  the  value  of  the  weighted  average  half  width.  A  typical 
value  of  the  opacity  grid  spacing  for  the  CH*  calculations  is  0.0005  cm*1  at  pressures  of  the  order 
of  1  mb.  The  normal  cutoff  point  in  the  line  profile  is  set  at  25  cm*1  from  line  center.  This  is  done 
both  for  reasons  of  economy  (beyond  25  cm1  a  given  line  contributes  little  opacity)  and  because 
deviations  from  Loren tz  behavior  can  become  appreciable.  In  addition  to  the  input  from  the 
SELECT  code.  The  values  of  temperature,  pressure,  and  nominal  mixing  ratios  (3Cs)  of  the 
absorbing  species  are  required  for  the  FASCODE  program.  These  values  are  obtained  from  a 
model  atmosphere  specification  which  will  be  discussed  later.  The  ABSORT  code  takes  the 
absorption  coefficient  files  generated  by  the  FASCODE  program  and  sorts  the  opacity  into  bins 
of  equal  logarithmic  width,  A  log  k(v),  to  produce  a  distribution, /ft),  based  on  the  relative 
probability  of  occurrence  within  the  wave  number  interval  (proportional  to  the  number  of  entries 
in  each  bin).  The  cumulative  distribution  function,  g(k)  (Equation  1),  is  obtained  by  numerical 
integration  of  the/?*)  function.  The  k-distributiou  is  obtained  by  a  tevetse  interpolation  of  the 
g(k)  relation  using  a  spline  function.  For  the  calculations  in  this  paper  a  201  bin  model  was  used. 
FP  suggests  a  number  greater  than  100  bins  for  normal  atmospheric  problems.  The  inputs  to  the 
ABSORT  code  ate  the  maximum  and  minimum  values  of  the  opacity  in  each  atmospheric  layer 
(determines  the  range  of  k’s  coveted  by  the  bins),  the  number  of  sorting  bins  (determines  the 
width  of  each  bin),  the  choice  Of  linear  or  spline  interpolation,  and  the  tine  by  tine  opacity  results 
from  FASCODE  for  each  layer.  The  output  from  the  ABSORT  code  is  the  201  point  k(g)  relation 
for  each  layer. 


; 


In  order  to  make  a  direct  comparison  with  the  LO  calculations,  the  wavelength  range 
chosen  for  the  CH4  and  NjO  calculations  in  this  paper  is  1100  to  1340  cm1  or  7.5  to  9.1  microns. 
This  wave  number  interval  was  divided  into  subintervals  having  a  20  cm1  width.  For  this 
subinterval  width  the  variation  of  the  value  of  the  Planck  function  is  less  than  10  percent  across 
the  band  at  the  maximum  and  minimum  atmospheric  temperatures,  Le.,  the  use  of  the  avenge 
value  of  the  Planck  function  over  the  wave  number  interval  is  accurate  to  the  order  of  about  5 


percent  of  the  actual  value.  The  use  of  wider  bands  will  require  averaging  opacities  weighted  by 
the  Planck  function 


The  model  atmosphere  chosen  for  the  calculations  is  the  McClatchey  mid-latitude,  summer, 
pressure-temperature  distribution  (McClatchey  a  al„  1972).  The  altitude  resolution  was  chosen 
as  1  km  between  0  and  20  km  and  2  km  between  20  and  62  km.  Mixing  ratios  of  CH4  and  tyO 
were  chosen  to  be  1.75  and  0.3  ppm  respectively.  Both  the  LO  calculations  and  those  of  this 
paper  assume  that  the  mixing  ratios  are  constant  with  altitude,  although  G1  indicates  that  the 
(nixing  ratios  decrease  significantly  above  -25  km  altitude.  This  point  will  be  discussed  later  in 
terms  of  the  stratospheric  cooling  rates. 

Figures  2  and  3  show  the /ft),  g(k),  and  k(g)  functions  for  CH4  in  the  wave  number  band 
1300-1320  cm*1  at  the  temperatures  and  pressures  representative  of  the  surface  layer. 

T  *  291.7  K  and  945.7  mb  (Figure  2),  and  a  layer  at  -50  km,  T  «=  272.2  K  and  P  =  1.48  mb 
(Figure  3).  It  is  interesting  to  note  that  the  fc-distributions  exhibit  very  different  behavior  at  g 
values  between  -0.9  and  1.0.  For  the  high  pressure  surface  layer  the  k  variation  is  less  than  an 
order  of  magnitude,  while  for  the  low  pressure  high  altitude  layer  the  k  variation  is  approximately 
four  orders  of  magnitude.  Figure  4  shows  die  ^-distributions  for  five  representative  atmospheric 
layers,  P  *  1.48, 18.05, 102.44, 542.36,  and  945.28  mb.  These  curves  are  constructed  from  the 
201  k(g)  points  produced  by  the  ABSORT  code.  For  pressures  greater  than  -100  mb  (altitudes 
below  -18  km)  the  curves  are  reasonably  well  behaved.  At  pressures  below  -100  mb  the  curves 
have  lane  opacity  excursions  at  g  value  greater  than  0.9.  These  excursions  can  be  up  to  four 
otden  of  magnitude  at  the  lowest  pressures.  This  kind  of  behavior  at  low  pressures  is  thought  to 
be  due  to  the  absence  of  pressure  brotdening  effects  on  the  absorption  lines  in  the  wave  number 
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band;  Le.  the  lines  ate  dominated  by  doppler  broadening  near  line  center.  Inspection  of  Figures  2 
and  3  shows  that  the  existence  of  narrow,  high  opacity  line  axes  in  the  wave  number  band 
causes  the  distribution  function  flk)  to  have  small  but  finite  values  at  high  opacities,  which 
results  in  a  sharply  peaked  k(g)  function  at  high  g  values.  These  variations  in  the  t -distributions 
require  a  careful  numerical  integration  strategy  in  the  transmission  expression,  Equation  3b,  in 
oner  to  accurately  reproduce  the  k(g)  functions.  The  integral  over  g  in  Equation  5b  can  be 
approximated  by  a  summation. 


N 

0  =  £fliexP 

«i  L  j 


Si  • 


(6) 


where  the  a,'  s  represent  the  integration  weights.  The  ^-distributions  for  the  lower  pressures 
require  a  larger  number  of  k(g)  points  in  Equation  6  to  accurately  reproduce  die  transmission 
integral.  The  integration  strategy  which  was  adopted  after  test  calculations  was  a  43  point, 
variable  spaced  trapezoidal  model.  For  g  values  between  0.0  and  0.6,  a  Ag  spacing  of  0.1  was 
chosen.  For  g  values  between  0.6  and  0.9,  a  Ag  spacing  of  0.02  was  chosen.  For  g  values 
between  0.9  and  1.0,  a  Ag  spacing  of  0.003  was  chosen  (this  is  the  spacing  for  the  201  point 
ABSORT  code  output  data). 

III.  Radiative  Transfer  Model 


In  order  to  determine  the  atmospheric  heating  cates  and  upward  and  downward  fluxes,  the 
radiative  transfer  equations  must  be  solved  for  die  specific  monochromatic  intensities, 
I„(Tv,/t),in  the  atmosphere.  The  classical  radiative  transfer  equations  for  a  plane  parallel 
atmosphere  are  (Liou,  1980); 


dxv 


U*i..m)-Bu  . 


dxv  =  ~kv  pdz  , 


(7a) 

(7b) 


H  ncosO  . 


(7c) 


ru  is  the  monochromatic  optical  depth  at  frequency  u(  Tv*0  at  the  top  of  the  atmosphere),  Bv  is 
the  Planck  function,  p  is  the  mass  density,  z  is  the  altitude,  and  0  is  the  zenith  angle.  The  formal 
solution  of  Equation  7a  is 


(8) 


For  downward  radiances,  -l£/i<0,  I,(0,-m)  =  0,  and  Equation  8  becomes 

(p>o)  . 


For  upward  radiances,  tu,  =  x„  (surface),  I v(rU(,p)«Bv  (ground),  l£fi<0,  and  Equation  8 
becomes 
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Integration  of  Equations  9  and  10  by  pans  gives 


for  the  downward  radiance  and 


lUT«.^)  =  Sv(fu)+[flu(8round)-fiu(fu,)V 

-|7 *•-’>  ^idr;  . 


for  the  upward  radiance.  Note  that  a  distinction  has  been  made  between  the  ground  temperature 
and  the  air  temperature  of  the  zero  altitude  point  The  flux,  F&u,  in  a  particular  frequency  band 
Au  is  given  by  the  equation 


Fl.  =  j  2xjl?(r¥,p)fu(fi. 


This  can  be  approximated  by  the  expression* 

Fav  -  j  I*.  l4(*v.  Mih  dv  , 

av  ‘*i  (14) 

where  the  a,’s  are  the  weights  for  the  numerical  integration  scheme.  Under  the  usual  assumption 
of  the  constancy  of  the  Planck  function  across  the  frequency  band  Av,  the  flux  equation  becomes 

Fa»“S>  /*( Ta,.A<,>i 

<= i  (15) 

where  TAv  is  the  transmission  function  for  the  frequency  interval  Av,  calculated  between  the 
appropriate  levels  in  the  atmosphere,  and  the  bi’s  are  the  integration  weights. 

The  /*  expressions  are; 

f  —  T"au(*'  zm*x’  Hi )"*■  j^4»  dB^v(z  ), 
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{*  =  SA¥(z)-[B4v(grourul)-fl4,(0)]TA,(z.o.^)-jTA,(z.2'>p1)^.(^). 

0  (17) 


and 


«*(*)- J*M*- 

A.  (18) 

This  set  of  equations  plus  the  transmission  equations  between  appropriate  levels  (cf.  Equation 
(6)),  are  solved  for  ten  evenly  spaced  p  values  between  0  and  1 .0  for  every  level  in  the 
atmosphere.  Once  the  upward  and  downward  fluxes  are  obtained  at  each  level,  the  heating  rates 
for  each  atmospheric  layer  can  be  determined.  For  a  particular  layer  the  heating  rate  is 

dT  AF 


where  AF  is  the  net  flux  into  the  layer, 

AF  =  F ufi  +  F  dnm  -  F -  F  DNi, 


(20) 


cf  is  the  specific  heat,  and  Am  is  the  mass  per  unit  area  in  the  layer.  The  subscripts  i  + 1  and  i 
denote,  respectively,  the  upper  and  lower  boundaries  of  the  layer.  The  hydrostatic  equilibrium 
condition  Am ■ pAz  =  -A Pig,  can  be  substituted  into  Equation  (19)  to  obtain  the  expression 


—(’Clday)  =  ~8-4672  ^  M 
dx  V  y>  AP  (mb)  ’ 

where  AP  is  the  pressure  difference  -P;)  across  the  layer. 

IV.  Results  and  Discussion 


(21) 


Upward  and  downward  fluxes  and  heating  rates  have  been  calculated  for  CH«  and  N?0 
using  die  correlated-*  transmission  model  for  the  following  problem  specifications: 

1.  Mcdatchey  mid  latitude,  summer,  model  atmosphere  temperature-pressure  distribution. 

2.  Altitude  resolution; 

a.  1  km,  0-20  km  altitude, 

b.  2  km,  20-60  km  altitude. 

3.  A  ground  temperature  of  294  K. 

4.  A  CHU  mixing  ratio  of  1.75  ppm,  constant  with  altitude. 

3.  An  N]0  mixing  ratio  of  0.31  ppm,  constant  with  altitude. 

6.  A  wavenumber  range  of  1100  to  1340  cm-*  (-7.5-9  microns),  in  20  cm*1  bands. 

Figures  3  (CH*)  and  6  (N2O)  show  die  downward  flux  (a),  the  heating  rate  (b),  the  upward 
flux  (c),  and  the  net  upward  flux  (d),  as  a  function  of  altitude  for  the  wavenumber  range 
1 100-1340  cm*1.  The  net  upward  flux  is  defined  as  the  upward  flux  at  each  atmospheric  level 
minus  the  upward  flux  due  to  the  ground  at  the  bottom  of  the  atmosphere.  The  net  upward  flux  is 
illustrated  for  the  purpose  of  direct  comparison  to  the  LO  results.  The  line  by  line  calculations  of 
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LO  for  the  downward  flux,  net  upward  flux,  and  heating  rates  are  shown  as  dotted  lines  in 
Figures  5a,  b,  d  and  6a,  b,  d.  The  wave  number  spacing  for  the  line  by  tine  calculations  was 
0.005  cm-1  for  altitudes  below  IS  km  and  0.0005 cm-1  for  altitudes  above  IS  km,  yielding 
between  4000  and  40000  points  within  the  20  cm*1  wave  number  band.  These  numbers  are  to  be 
compared  with  the  43  k(g)  points  for  the  same  20  cm*1  wave  number  band  for  the  k  -distribution 
transmission  model  The  results  for  CH4  (Figure  5)  show  that  fluxes  agree  10  -8  percent  with  the 
line  by  line  calculations.  The  heating  rate  calculations  show  good  agreement  up  to  an  altitude  of 
-45  km  with  the  line  by  line  results.  The  quality  of  the  agreement  below  45  km  is  of  the  same 
order  as  that  given  by  the  LO  results.  Above  45  km  the  heating  rate  shows  a  large  negative 
excursion  which  is  correlated  with  a  temperature  maximum  at  47  km.  In  the  region  between  45 
and  60  km  the  * -distribution  model  produces  heating  rates  which  are  greater  than  the  line  by  line 
results  by  -4  percent  Test  calculations  seem  to  indicate  that  the  agreement  of  die  heating  rates  in 
these  altitude  regimes  depends  on  how  accurately  the  high  opacity  end  of  the  tow  pressure  k- 
distributions  can  be  resolved. 

The  results  for  NjO  (Figure  6)  exhibit  similar  behavior  to  the  CHU  results.  The  fluxes  agree 
with  the  line  by  line  values  to  better  than  -15  percent.  The  heating  rates  below  -45  km  show 
good  agreement  with  the  line  by  line  values  and  is  of  the  same  quality  as  die  CH4  results.  Above 
45  km  there  is  a  disagreement  in  the  heating  rases  which  amounts  to  -5  percent 

As  a  test  of  the  sensitivity  of  the  method  to  the  size  of  the  wavenumber  subintervals,  n  set 
of  heating  rate  calculations  were  made  using  a  subinterval  size  of  40  cm4.  The  headng  rate 
variation  as  a  function  of  altitude  at  the  40  cm*1  interval  size  for  CH«  and  N2O  is  shown  in 
Figures  7a  and  7b  respectively.  The  accuracy  of  die  40  car1  interval  calculation  in  the  45 10 
60  km  altitude  range  decreases  compared  to  20  cm*1  interval  size  calculation.  For  CH4  the  error 
in  the  heating  rate  goes  from  -4  percent  to  ~12  percent  as  the  interval  size  goes  from  20  cor1  to 
40  cm-1.  For  NjO  the  enor  in  the  heating  rate  goes  from  -5  percent  to  -42  percent  as  the  interval 
size  increases  from  20  cm*1  to  40  cnr*.  The  decrease  in  the  hearing  rate  error  is  due  to  two 
factors,  a  smaller  variation  of  the  Planck  function  across  the  interval  and  increase  in  resolutioa  of 
the  ^-distribution.  The  allocation  of  the  error  difference  between  these  two  factors  depends  on  the 
spectral  characteristics  of  the  gases  being  considered. 

The  CH4  and  N2O  heating  rate  calculations  of  LO,  using  Jt -distributions  based  on  the 
Malkmus  Band  Model  to  calculate  the  transmissions,  agree  with  the  line  by  line  results  to  within 
approximately  5  percent  on  the  45  to  60  km  altitude  region  (since  this  is  the  alrimde  region  when 
doppler  line  cores  become  important).  The  agreement  at  the  lower  altitudes  is  much  better.  In  the 
0-60  km  altitude  region,  the  FP  Jb-distribution  model  used  in  this  paper  and  the  Malkmus  k- 
distriburion  model  used  by  LO,  produce  approximately  die  same  results.  The  calculations  in  this 
paper  have  been  done  using  a  wave  number  band  spacing  of  20  cm'1  and  43  k  intervals  for  the  k- 
distribudon  in  order  to  resolve  the  sharply  peaked  k  distributions  in  the  stratosphere.  These 
numbers  are  to  be  compared  with  the  LO  band  spacing  of  10  cm*1  and  36  k  intervals  for  the  k- 
distribution. 

The  results  of  the  CH«  and  N2O  calculations  indicate  that  the  correUted-k  method  of 
calculating  transmissions  appears  to  be  a  computationally  efficient  method  which  yields  good 
accuracies  when  compared  with  line  by  line  calculations.  Future  work  will  address  the  overlap 
capabilities  of  the  £ -distributions  for  different  species  as  well  as  the  inclusion  of  more  molecules 
in  order  to  produce  a  complete  longwave  tadiaoon  transport  model  Additional  efforts  will 
include  numerical  improvements  to  optimize  quadratures,  allow  wider  frequency  bands,  the 
implementation  of  a  new  code  to  calculate  the  finely  gridded  set  of  monochromatic  absorption 
coefficients  used  by  the  ABSORT  code,  and  the  addition  of  clouds  and  aerosols. 
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Figure  Captions 

Figure  1.  A  block  diagram  of  the  calculation  of  the  upward  and  downward  fluxes  and  the 

healing  rates.  The  SELECT  Code  accesses  the  HITRAN-91  data  base  and  outputs  the 
appropriate  line  transition  data  for  the  selected  molecular  species.  The  FASCODE 
Code  takes  the  output  of  the  SELECT  Code  plus  the  temperature-pressure  relation 
and  mixing  ratio  values  from  the  model  atmosphere  and  produces  a  line  by  line 
opacity  vs.  wave  number  data  base-far  the  selected  wave  number  band.  The  ABSORT 
Code  uses  the  line  by  line  data  base  output  from  the  FASCODE  Code  to  calculate  the 
fik),  g(k),  and  k(g)  relations  for  each  atmospheric  point  The  CORK1  Code  uses  the 
k(g)  relations  calculated  in  the  ABSORT  Code  and  the  model  atmosphere  to  calculate 
fluxes  and  heating  rates  at  each  atmospheric  level. 

Figure  2.  The  CH4  opacity  distribution  function,  fik)  (Figure  2a),  cumulative  probability,  gik ) 
(Figure  2b),  and  4 -distribution,  k(g)  (Figure  2c),  for  the  wave  number  band  1300- 
1320  cm-1,  at  T  =  272.23  K  and  P  =  1.48  mb. 

Figure  3.  The  CH4  opacity  distribution  function,  fik)  (Figure  3a),  cumulative  probability.  g(k) 
(Figure  3b),  and  * -distribution,  k(g)  (Figure  3c),  for  the  wave  number  band  1300- 
1320  cm1,  at  T  -  291.66  K  and  P  =  945.2  mb. 

Figure  4.  The  CH4  ^-distributions  for  the  wave  number  band  1300-1320  cm1,  at  five 
atmospheric  pressure  levels,  P  =  1.48, 18.05, 102.44, 542.36,  and  945.28  mb. 

Figure  5.  The  altitude  variation  of  the  radiative  fluxes  and  heating  rate  for  CH4  in  the  wave 
number  band  1100-1340  cm*1  using  a  29  cnr1  subinterval  size.  Figure  5a  shows  the 
altitude  variation  of  the  downward  flux  (W/m2).  Figure  5b  shows  the  altitude 
variation  of  the  heating  rate  (deg  C/day).  Figure  5c  shows  the  altitude  variation  of  the 
upward  flux  (W/m2).  Figure  5d  shows  the  altitude  variation  of  the  net  upward  flux 
(W/m2).  The  net  upward  flux  is  defined  as  the  upward  flux  at  a  particular  atmospheric 
level  minus  the  upward  flux  due  to  the  ground  at  the  bottom  of  the  atmosphere.  The 
dotted  lines  shown  in  Figures  5a,  Sb,  and  Sc  are  the  line  by  line  values  of  LO. 

Figure  6.  The  altitude  variation  of  the  radiative  fluxes  and  heating  rate  for  N2O  in  the  wave 
number  band  1 100-1340  cm*1  using  a  29  cnr1  subinterval  size.  Figure  6a  shows  the 
altitude  variation  of  the  downward  flux  (W/m2).  Figure  6b  shows  the  altitude 
variation  of  the  heating  tale  (deg  C/day).  Figure  6c  shows  the  altitude  variation  of  the 
upward  flux  (W/m2).  Figure  6d  shows  the  altitude  variation  of  the  net  upward  flux 
(W/m2). The  net  upward  flux  is  defined  as  the  upward  flux  at  a  particular  atmospheric 
level  minus  the  upward  flux  due  to  the  ground  at  the  bottom  of  the  atmosphere.  The 
dotted  lines  shown  in  Figures  6a,  6b,  ami  6c  are  the  line  by  line  values  of  LO. 

Figure  7.  The  altitude  variation  of  die  heating  rate  (deg  C/day)  in  the  wave  number  band  1 100- 
1340  cm*1,  using  a  40  cnr1  subinterval  size.  Figure  7a  shows  the  CH<  heating  rate 
variation  and  Figure  7b  shows  the  tyO  heating  rate  variation. 
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LINE  BY  LINE  CALCULATION  OF  ATMOSPHERIC  FLUXES  AND  COOLING 
RATES:  APPLICATION  TO  WATER  VAPOR,  OZONE,  CARBON  DIOXIDE 

AND  THE  FLUOROCARBONS 

SA.  Clough,  MJ.  Iacono,  J.-L.  Moncet 
Atmospheric  and  Environmental  Research,  Inc., 

840  Memorial  Drive,  Cambridge,  MA  02139 

Line  by  line  calculations  of  atmospheric  fluxes  and  heating  rates  for  die  long  wave  spectral 
region  have  been  performed  to  assess  methods  of  treating  the  overlap  of  caibon  dioxide  and  water 
vapor  and  to  study  approximate  methods  for  including  the  effects  of  heavy  molecules  including  the 
fluorocaibons.  In  developing  these  results  we  have  explicitly  considered  the  role  of  the  water 
vapor  continuum  [1]  and  the  effect  of  doubling  carbon  dioxide.  The  role  of  die  water  vapor 
continuum  is  explicitly  considered.  Comparisons  with  results  from  ICRCCM  will  be  presented.  A 
vectorized  version  of  FASGODE  with  the  capability  of  obtaining  fluxes  has  been  utilized  far  these 
calculations.  The  strong  conelation  of  die  spectral  cooling  rates  and  the  outgoing  spectral  radiance 
is  emphasized  as  well  as  the  importance  of  the  spectral  region  from  0  -  600  cm'1. 

1 .  Clough,  S.A.,  F.X.  Kneizys  and  R.W.  Davies,  1989:  line  Shape  and 
the  Water  Vapor  Continuum,  Atmospheric  Research,  22. 229-241. 
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BACKSCAT  LIDAR  SIMULATION:  VERSION  3.0 

J.R.  Hummel,  D.R.  Longtin,  N.L.  DePiero,  RJ.  Gras  so 
SPARTA,  Inc.,  24  Hartwell  Avenue,  Lexington,  MA  02173 

The  Geophysics  Directorate  of  Phillips  Laboratory  is  developing  a  number  of  lidar  systems 
for  use  in  probing  die  atmosphere.  These  systems  include  backscatter  lidars  to  study  atmospheric 
aerosols,  Doppler  lidar  systems  to  measure  wind  fields,  and  Raman  lidars  to  study  the  distributions 
of  different  molecular  species.  To  aid  in  the  design  and  use  of  such  lidar  systems,  SPARTA  has 
developed  a  lidar  simulation  program,  BACKSCAT.  Originally  developed  to  include  only  die 
badacattered  return  from  aerosols,  the  simulation  package  has  evolved  to  also  include  Raman 
scattering  processes.  This  paper  describes  a  new  version  of  the  simulation  system,  BACKSCAT 
Version  3.0. 

BACKSCAT  Version  3.0  includes  two  significant  improvements.  The  first  is  die  inclusion 
of  user-defined  aerosol  layers  and  die  second  is  die  consideration  of  Raman  scattering  processes. 
(The  description  of  the  Raman  scattering  options  within  BACKSCAT  Version  3.0  will  be 
described  in  a  companion  paper.) 

In  BACKSCAT  Version  3.0,  a  user-defined  aerosol  layer  is  defined  by  a  number  density 
profile,  a  sire  distribution  shape,  and  an  index  of  refraction.  Aerosol  attenuation  properties  are 
computed  using  an  efficient  Mie  scattering  program  that  is  coupled  to  the  BACKSCAT  simulation 
system.  Users  can  select  from  a  library  of  aerosol  indices  of  refraction  for  common  aerosols,  or 
they  can  input  specific  values. 

Research  Supported  by  Phillips  Laboratory,  Geophysics  Directorate,  Contract  F19628-91-C-0093. 


288 


L  A 


OVERVIEW  OF  BACKSCAT  VERSION  3.0 
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From  a  Field  Program 


EVOLUTION  OF  THE  BACKSCAT 


-  User-Defined  Aerosols 

-  System  Efficiency  Considered 

-  Raman  Lidars  Simulated 


OVERVIEW  OF  BACKSCAT  VERSION  3.0 


Molecular  Absorption 
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•  System  Efficiency  Added 
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Full  Single-to-Noise  Treatment 
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NEW  FEATURES 
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-  Particle  Complex  Index  of  Refraction 
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Smoke 

Maritime  Aerosols 
Stratospheric  Aerosols 
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01  0.001  , 
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-  Additions  of  a  Raman  Lidar  Option 
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•  Include  the  Capability  to  Invert  Solutions  From 
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DEVELOPMENT  OF  A  RAMAN  LIDAR  SIMULATION  OPTION  WITHIN 

BACKSCAT  VERSION  3.0 

R.J.  Gras  so,  J.R.  Hummel 

SPARTA,  Inc.,  24  Hartwell  Avenue,  Lexington,  MA  02173 

Raman  lidar  is  a  useful  and  powerful  tool  for  remote  probing  of  the  a&nosphere.  With 
Raman  lidars,  one  can  accurately  determine  die  identification  and  concentration  of  a  particular 
molecular  specie  present  in  die  atmosphere.  Results  are  presented  from  an  effort  to  develop  a 
simulation  capability  of  Raman  lidar  systems  for  the  remote  detection  of  atmospheric  gases  and/or 
air  polluting  hydrocarbons.  Our  model,  which  integrates  remote  Raman  spectroscopy  with 
SPARTA's  BACKSCAT  Version  3.0  atmospheric  lidar  simulation  package,  permits  accurate 
determination  of  die  performance  of  a  Raman  lidar  system. 

The  accuracy  of  the  model  results  from  the  accurate  calculation,  at  any  given  excitation 
wavelength,  of  foe  differential  scattering  cross  section  for  the  molecular  specie  under  investigation. 
We  show  excellent  correlation  of  our  calculated  cross  section  data  with  experimental  data  from  the 
published  literature.  In  addition,  the  use  of  BACKSCAT  Version  3.0  package,  which  provides  a 
user  friendly  environment  to  define  die  operating  conditions,  provides  an  accurate  calculation  of  the 
atmospheric  extinction  at  both  the  excitation  and  Raman  shifted  wavelengths.  The  Raman  option 
within  BACKSCAT  Version  3.0  can  be  used  to  accurately  predict  the  performance  of  a  Raman 
lidar  system,  the  concentration  and  identification  of  a  specie  in  the  atmosphere,  or  the  feasibility  of 
making  Raman  measurements. 

Research  Supported  by  Phillips  Laboratory,  Geophysics  Directorate,  Contract  F19628-91-C-0093 
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DEVELOPMENT  OF  A  RAMAN  LIDAR 
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Annual  Review  Conference  on  Atmospheric  Transmission  Models 

Hanscom  AFB,  MA  01731-5000 

Funded  By  Air  Force  Contract  FI  9628-91 -C-0093 


OUTLINE  OF  BRIEFING 
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Implementing  a  Raman  Simulation  in  the  BACKSCAT  Version  3.0  System 


BACKGROUND  AND 
PURPOSE  OF  THE  STUDY 
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Implement  a  Raman  Lidar  Simulation  Capability  Within  the  BACKS 
Version  3.0  Framework 


OVERVIEW  OF  THE  RAMAN  EFFECT 
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IMPLEMENTING  A  RAMAN  SIMULATION 
IN  THE  BACKSCAT  FRAMEWORK 
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Solid  Angle  and  Ni(R)  is  the  Concentration  of  the  Species  Being  Probed 
T(X0,  R)  and  r(Ar,  R)  are  the  One-Way  Total  Atmospheric  Transmissions 
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Most  Cross  Section  Data  Measu 
Lasers  at  488  nm  and  514.5  nm 
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Wavelengths  Shorter  Than  337  nm 


SUMMARY  OF  RAMAN  CROSS  SECTIONS  FOR 

NITROGEN 
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RELATIVE  CROSS  SECTIONS  FOR 
C02,  H20  VAPOR,  03,  &  02 


RAMAN  CROSS  SECTION  MODEL 
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RAMAN  CROSS  SECTION  MODEL 

(Cont.) 
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Tensor 

h  is  Planck’s  Constant,  k  is  Boltzmann’s  Constant,  c  is  the  Velocity  of  Light, 
and  T  is  the  Temperature 
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is  Independent  of  the  Scattering  Geometry 


RAMAN  CROSS  SECTION  MODEL 
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For  Molecules  Where  No  Previous  Raman  Data  Exists,  ^(45 
Would  Have  to  be  Calculated 


250  300  350  400  450  500  550  600  650  700 

Wavelength  (nm) 


RESULTS  FROM  RAMAN  CROSS  SECTION 
MODEL  -  CARBON  DIOXIDE 


Wavelength 


SUMMARY  AND  RECOMMENDATIONS 
FOR  FUTURE  WORK 
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User  Profiles 


SUMMARY  AND  RECOMMENDATIONS 

(Cont.) 
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L  A 


CALCULATION  OF  THE  ANGULAR  RADIANCE  DISTRIBUTION  FOR  A 
COUPLED  SYSTEM  OF  ATMOSPHERE  AND  CANOPY  MEDIA  USING  AN 
IMPROVED  GAUSS-SEIDEL  ALGORITHM 


S.  Liang,  A.H.  Strahler 


onwealth  Avenue,  Boston,  MA  02215 


The  radiative  transfer  equations  of  coupled  system  of  atmosphere  and  canopy  media  are 
solved  numerically  by  an  improved  Gauss-Seidel  iteration  algorithm.  Hie  radiation  field  is 
decomposed  into  three  components:  uncollided  sunlight,  single  scanning,  and  multiple  scattering 
radiance  far  which  the  corresponding  equations  and  boundary  conditions  are  set  up  and  their 
analytical  or  itenuional  solutions  are  explicitly  derived.  The  hotspot  effect  of  the  canopy  is 
accomplished  by  means  of  the  modification  of  extinction  coefficients  of  upiwud  single  scattering 
radiation  and  uncoBided  sunlight  To  reduce  the  computation  for  the  case  of  deo>  optical 
thickness,  an  improved  iteration  formula  is  derived  to  speed  convergence.  The  formulation 
presented  in  this  paper  is  well  suited  to  analyze  the  relative  magnitude  of  multiple  scattering 
radiance  and  single  scattering  radiance  in  both  the  visible  and  near  infrared  regions. 


r 


325 


r  i 


CALCULATION  OF  THE  ANGULAR  RADIANCE  DISTRIBUTION 
FOR  A  COUPLED  ATMOSPHERE  AND  CANOPY 


Shunlin  Liang  &  Alan  H.  Strahler 
Center  for  Remote  Sensing,  Boston  University 


1.  Introduction 


2.  Radiative  Transfer  Models 

*  atmospheric  radiative  transfer  equation 

*  canopy  radiative  transfer  equation 

*  canopy  leaf  model 


3.  An  Improved  Gauss-Seidel  Numerical  Algorithm 


4.  Model  Validation  and  Data  Analysis 

5.  Conclusions 
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Introduction 


*.  Limitations  of  the  conventional  remote  sensing 

based  on  nadir  observation 

\ 


*.  Spaced 

sible  in  EOS  era 


multiple  angle  observations  are  pos- 


*.  Individual  modeling  efforts  for  the  atmosphere 
and  canopy 
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L  A 


Atmospheric  radiative  transfer  equation 

=n%,a)- ~falj>(ar->Gmx.c*)dcr 

Canopy  radiative  transfer  equation 

-  + h  (t,  Q)  g  (fl)  /  (t,  q> = i  J^r(n,->C2)  /  (t,  nvn' 

subject  to  corresponding  boundary  conditions. 
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The  function  G(fl)  is  the  mean  projection  of  a  unit  foliage  area 
in  the  direction  fi,  i.  e. 


G(Q)  = 


The  correction  function  for.  hotspot  effect  is  given  by 


A( x,Q)=l 


A(Qo.Q)x 
— JR — 


It  can  be  observed  that  in  the  case  of  backscattering  (i.  e., 
Q=-fto),  A(Qo,0)  =  0  and  h(x,ii)= 0.  The  absence  of  extinction 
coefficient  results  in  the  local  maximum  of  reflectance  which  is 
widely  termed  the  hotspot  peak. 
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L  A 


The  radiative  transfer  equations  and  their  boundary  conditions 
are  given  by 

-ni£^LQI+/*(t.Q)= p(Q'->ny"(t,«vor+e,(t.Q) 

- +  g (tiy * (t, Q) = -i- j^r(0'-*Qy* (x, ova' + G2(*. O)  T- <Tit* 

/*(0,Q)  =  0  |x<o 

/*f(x/,Q)=  f  |i  >  0 

“■  l*  ,  ' 


Q  i(t.Q) = ^-j^(a'->nyO(x,QVQ' + £  J4  *(x,QW 1 5  xa 

q2(x,o) = -^4^  °(x.  over + ^r(Q,-»a)/  *(x.  ayt  o'  T»  <TSx' 

> 

For  simplicity  of  discussions,  assume  that  we  are  faced  with 
such  a  general  equation 

-pPJU&Sb  +  /(Q)/«(x,Q)  =7(x,Q) 


where 


f(0) 


_fl  x^xtt 
-|C(£1)  xa  <x£x, 


and  the  source  function  is 


^-J^p (Q’->n)[iM(T, £T) + / «(x. n' )]</n'  +  £j^p<a^y<to.aw  T^ttf 
lj  r(O'-»QK/A#(x.o')+/,(x,n,)wn,+i-4^r(Q'-^ii)/0(x,£i')dn'  %a  <t5 


Gauss-Seidel  algorithm 


From  die  spatial  interval  Herman  and  Biowning  assume 

die  source  function  J(x,  Q)  may  be  tatcen  to  be  independent  of  t'  and  equal 
to  its  value  at  die  midpoint  of  die  interval,  toting  /(t'.Q)  =J (tj+t,Q)  for 
all  t*  in  this  interval. 

The  improved  version  is  base  on  a  linear  variation  of  die  source  func¬ 
tion.  For  three  arbitrary  layers  t,-  £  t'  £  t(+2,  the  linear  relation  can  be  given 

by: 

=j(xi+2,Q)  +  (t, _ w  ]l>0 

For  the  downward  radiance, 

/*(t<+2,Q)  =  /^,12)esp(-2At')  + 

-:^^-[l-,/Ax'  +  (1+At,)exp(-2At')]  p<0 

The  same  procedures  will  yield  the  formula  for  the  upwelling  direction 
/^(x,  ,Q)  =  /^(xl+2,«)exp(-2AT/)  + 

T^T'fl  ~  +d  +A^)exp(-2AT')J  p>0. 
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Tfcbte  Validation  of  Gntm-Seiiei  Aijoritha  Ifoiaf  DbcrMe-OtdiaMe  Alcoruhta 


viswiag  aagta 

t=0.1 

T»1A 

t 

=3j0 

|iscom 

Qmm-SMd 

DtaaetoOidiBMB 

Om -Seidel 

Dta*°*“e 

0067 

043419 

043408 

059566 

059400 

062154 

0383 

035950 

03S95S 

047818 

047954 

052799 

053298 

0574 

Bi  lit 

032146 

032526 

032689  | 

0J9136 

059783 

0.840 

031234 

024526 

034652 

03926S 

039904 

0987 

030962 

030963 

031330 

031441 

034804 

034803 

incident  net  flax  /'o=K,  *otl  reflectance  /?,  =0.3,  Ho =0.70292,  $= 60.  00=0.96.$  =0.65. 


32-cnen  far  dbci«)o<m&MiB  dgarifl^  30  by  20  far  Oaan^ddei  dfaiidn. 
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Conclusion 

*.  The  improved  Gauss-Seidel  algorithm  can  be 
effectively  used  for  the  coupled  medium 

*.  We  could  not  accurately  retrieve  canopy  informa¬ 
tion  in  visible  region  without  exact  knowledge  of 
atmospheric  characteristics 

*.  The  near-IR  region  is  of  great  significance  for  the 
remote  sensing  of  vegetation 


L 


337 


A  PHYSICALLY  REASONABLE  ANALYTIC  EXPRESSION  FOR  THE  SINGLE 

SCATTERING  PHASE  FUNCTION 

W.M.  Comette,  J.G.  Shanks 

Photon  Research  Associates,  Inc.,  9393  Towne  Centre  Drive,  Suite  2000,  San  Diego,  CA  92121 

An  analytic  phase  function  that  reduces  die  Rayleigh  phase  function  for  die  scanning  of 
unpolarized  light  is  presented  and  compared  to  the  traditional  Henyey-Greenstein  phase  function. 
Comparisons  between  the  proposed  phase  function  and  the  phase  function  for  three  of 
Deirmendjian's  polydispersions  are  shown,  and  applications  to  radiative  transfer  are  demonstrated. 
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Dfego,  CaRfomte  *2121 
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NEW  PHASE  FUNCTION 
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ANALYSIS  OF  SOLAR  TRANSMISSION  DATA  OVER  THE  SOUTH  ATLANTIC 

OCEAN  DURING  SABLE  89 

D.  Longtin,  J.  Hummel 

SPARTA,  Inc.,  24  Hartwell  Avenue,  Lexington,  MA  02173 

G.G.  Koenig 

PL/GDOPA,  Hanscom  AFB,  MA  01731-5000 

This  paper  discusses  the  analysis  of  solar  transmissometer  data  that  was  taken  during  the 
South  Atlantic  Backscatier  Iidar  Experiment  (SABLE).  There  were  six  days  when  the  instrument 
was  operated:  25  June,  28  June,  1  July,  2  July,  6  July,  and  7  July  1989.  The  instrument  was 
loaded  on  Ascension  Island,  and  it  measured  solar  transmission  at  532  nm  versus  time  and 
transmission  across  the  solar  spectrum. 

Ia  the  analysis,  die  time  series  of  solar  transmissions  at  532  nm  are  used  to  study  the 
properties  of  cirrus  and  boundary  layer  stratocumulus  that  passed  in  front  of  the  sun  while  the 
instrument  was  working.  For  25  June  1989,  a  qualitative  assessment  shows  the  transmissions 
through  cirrus  to  be  between  0.4  and  0.8,  and  the  cirrus  thicknesses  are  estimated  to  be  between 
1.0  and  2.2  km.  These  estimates  arc  consistent  with  cirrus  observations  from  aircraft 

The  remainder  of  the  analysis  focuses  on  the  properties  of  stratocumulus  cloud  edges  and 
thin  spots.  It  is  shown  that  stratocumulus  have  very  similar  edge  characteristics  from  day-to-day: 
there  is  a  uniform  distribution  of  cloud  transmission  values  between  about  0.1  and  0.7,  and  the 
percentages  increase  significantly  as  cloud  transmissions  increase  from  0.7  to  0.9.  Also,  the 
number  of  occurrences  of  stratocumulus  thin  spots  and  edges  are  compared  against  total  number  of 
stratocumulus  occurrences.  Although  day-to-day  variations  exist,  the  comparisons  clearly  show 
that  stratocumulus  over  Ascension  Island  often  are  not  opaque,  and  they  frequently  transmit  partial 
solar  radiation.  Since  stratocumulus  often  cover  broad  regions  of  the  marine  environment,  these 
findings  could  be  of  interest  to  climate  modeling  studies  which  often  assume  stratocumulus  to  be 
fully  opaque. 
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-  Solar  Scans  Between  400  and  1200  nm  (Not  Used  in  This  Analysis) 
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